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Benchmark  datasets 

Validated  codes  for  subcomponent  applications 
CFD-derived  Advanced  Hardware  Concepts 
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- Increased  predicted  efficiency  over  conventional  designs  by  2 percent 


ANALYSIS  OF  A NEW  TURBINE  DESIGN 

Instant, moons  Mach  Numbor  Contours 
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Team  currently  exploring  mechanisms  to  reduce  losses  in  the  Gas  Generator  Oxidizer  Turbine 

— Minishrouds 

--  Blade  and  endwall  fences 

— Labyrinth  effects 
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Design,  build,  and  test  advanced  Gas  Generator  Oxidizer  Turbine  system  model 

- Design  verification 

— Code  validation 
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Verified  advanced  impeller  design  experimentally 


11  INCH  SSME  HPFTP  IMPELLER 

LASER  VFLOCIMETER  TEST  DATA 
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VT28  - Hydrodynamic  Design  of  Advanced  Pump  Components  TURBOMACHINERY 

Particle  Trace  of  the  Flow  Field  technology 

- Original  Design 
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Hub  Surface  Mid  Span  Tip  Surface 
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Performed  a first  order  impeller-diffuser  interaction  prediction 
- Used  existing  Caltech  data 
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Expansion/deflecting  nozzles 
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Innovative  design,  manufacturing,  and  management 
to  lower  cost 

Expect  NRA  to  be  released  4Q-93 
Henry  Stinson,  EP62 
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ARPA  looking  for  ’in-kind'  contributions 

Cannot  count  NASA  cooperation  as  part  of  'industry  matching  funds 
Dollars  contributed  by  NASA  subtract  from  ARPA  funding 
Use  of  NASA  facilities  and  manpower  does  not  affect  ARPA  funding 
(may  be  viewed  as  favorable) 
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- Low  Turbine  Temperatures 
LOW  COSTS 

- 50K  Upper  Stage  Engine  < $1M 

- Simple  Design  Supports  Minimum  Development  Program 


KEY  PERFORMANCE  PARAMETERS 
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- Physical  models  and  implementation 

- Benchmarking 

- Program-supporting  applications 

- Integration  of  CFD  into  design  organizations 
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unstructured  solvers 

adaption  techniques  and  error  estimates 

algorithm/ hardware  architecture  optimization 


National  Aeronautics  and  Space  Administration  ^ , 

George  C.  Marshall  Space  Flight  Center  OverView  Of  MSFC 

Structures  and  Dynamics  Laboratory  - • « . 

Fluid  Dynamics  Division,  CFD  Branch  LrD  Activities 
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Expanded  opportunities  to  integrate  CFD  into  advanced  space 
programs,  applications  for  dual-use  technologies. 


N95- 23616 
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3-D  CFD  ANALYSIS  OF  HYDROSTATIC  BEARINGS 

:)  i <3 

Shyi-Jang  Lin  and  Robert  I.  Hibbs  / - ' 

Rocketdyne  Division,  Rockwell  International 


ABSTRACT 

The  hydrostatic  bearing  promises  life  and  speed  characteristics  currently 
unachievable  with  rolling  element  bearings  alone.  In  order  to  achieve  the 
speed  and  life  requirements  of  the  next  generation  of  rocket  engines, 
turbopump  manufacturers  are  proposing  hydrostatic  bearing  to  be  used  in  place 
of,  or  in  series  with,  rolling  element  bearings. 

The  design  of  a hydrostatic  bearing  is  dependent  on  accurate  prediction  of  the 
pressure  in  the  bearing.  The  stiffness  and  damping  of  the  hydrostatic  bearing  is 
very  sensitive  to  the  bearing  recess  pressure  ratio.  In  the  conventional 
approach,  usually  ad  hoc  assumptions  were  made  in  determining  the  bearing 
pressure  of  this  approach  is  inherently  incorrect. 

In  the  present  paper,  a more  elaborate  approach  to  obtain  the  bearing  pressure 
is  used.  The  bearing  pressure  and  complete  flow  features  of  the  bearing  are 
directly  computed  by  solving  the  complete  3-D  Navier-Stokes  equation  . 

The  code  used  in  the  present  calculation  is  a modified  version  of  REACT3D 
code. 

Several  calculations  has  been  performed  for  the  hydrostatic  bearing  designed 
and  tested  at  Texas  A&M.  Good  agreement  has  been  obtained  between 
computed  and  test  results.  Detailed  flow  features  in  the  bearing  will  be  also 
described  and  discussed. 
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3-D  CFD  Analysis  of  Hydrostatic  Bearings 
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Benefits  of  Hydrostatic  Bearings  in  High  Power  Density 
Turbomachinery  Led  Rocketdyne  to  Pursue  Aggressive  IR&D 
Initiative  to  Improve  Analysis  Capability 


HPdTP  Pump  End  Bearing  Conversion 
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FUEL  PUMP  FOR  NLS 
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Rockwell  International 

Rockstdyne  Division 


Significantly  Improved  Life  & Power 
Level  Margin  With  Hydrostatic  Bearing 
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Rockwell  International 


Analysis  Method: 
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• Steady  Solution  with  Eccentric  bnan 

• Unsteady  Solution  with  Whirling  Shaft 


COMPARISON  OF  THEORY  AND  EXPERIMENT 

TEXAS  A&M  HYDROSTATIC  BEARING  TESTING 

Difference 

MIN  AVG  MAX 
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Rockwell  International 


3D  Mult:~Zcne Grid  off H$S  Bearing 
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Rockwell  International 

Rocketdyne  Division 
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Rockwell  International 

Rocketdyne  Division 
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Rockwell  International 

Rocketdyne  Division 


Code  Development  and  Verification 

TAMU  Data  Compared  to  CFD  Solution 
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Rockwell  International 

Rocketdyne  Division 


Pressure  Contour  at  Clearance  Radial  Plane 
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Rockwell  International 

Rockatdyne  Division 
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Finite-difference  Solutions  of  the  Alternate  Turbopump  Development  High- 
pressure  Oxidizer  Turbopump  Pump-end  Ball-bearing  Cavity  Flows 

by 


Theodore  G.  Benjamin,  Roberto  Garcia,  Paul  K.  McConnaughey, 
Ten-See  Wang  and  Bruce  T.  Vu 
Computational  Fluid  Dynamics  Branch 
George  C.  Marshall  Space  Flight  Center 
National  Aeronautics  and  Space  Administration 
Huntsville,  Alabama 
and 

Youssef  Dakhoul 
Sverdrup  Technologies 
Huntsville,  Alabama 


These  analyses  were  undertaken  to  aid  in  the  understandng  of  flow 
phenomena  in  the  Alternate  Turbopump  Development  (ATD)  High-pressure 
Oxidizer  Turbopump  (HPOTP)  Pump-end  ball  bearing  (PEBB)  cavities  and  their 
roles  in  turbopump  vibration  initiation  and  bearing  distress.  This  effort  was  being 
performed  to  provide  timely  support  to  the  program  in  a decision  as  to  whether  or 
not  the  program  should  be  continued. 

In  the  first  case,  it  was  determined  that  a change  in  bearing  throughflow  had 
no  significant  effect  on  axial  preload.  This  was  a follow-on  to  a previous  study 
which  had  resulted  in  a redesign  of  the  bearing  exit  cavity  which  virtually 
eliminated  bearing  axial  loading. 

In  the  second  case,  a three-dimensional  analysis  of  the  inner-race-guided 
cage  configuration  was  performed  so  as  to  determine  the  pressure  distribution  on 
the  outer  race  when  the  shaft  is  0.0002"  off-center.  The  results  indicate  that 
there  is  virtually  no  circumferential  pressure  difference  caused  by  the  offset  to 
contribute  to  bearing  tilt. 

In  the  third  case,  axisymmetric  analyses  were  performed  on  an  outer-race 
guided  cage  configuration  to  determine  the  magnitude  of  tangential  flow  entering 
the  bearing.  The  removed-shoulder  case  was  analyzed  as  was  the  static  diverter 
case.  A third  analysis  where  the  preload  spring  was  shielded  by  a sheet  of  metal 
for  the  baseline  case  was  also  performed.  It  was  determined  that  the  swirl 
entering  the  bearing  was  acceptable  and  the  project  decided  to  use  the  outer- 
race-guided  cage  configuration. 

In  the  fourth  case,  more  bearing  configurations  were  analyzed.  These 
analyses  included  thermal  modeling  so  as  to  determine  the  added  benefit  of 
injecting  colder  fluid  directly  onto  the  bearing  inner-race  contact  area.  The 
results  of  these  analyses  contributed  to  a programmatic  decision  to  include 
coolant  injection  in  the  design. 
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20  April  1993 


. «...  . . Computational  Fluid  Dynamics  Branch 

Finite-difference  Solutions  of  the  Alternate  K Rujd  Dynamics  Division 

Turbopump  Development  High-pressure  Oxidizer  Structures  and  Dynamics  Laboratory 

Turbopump  Pump-end  Ball-bearing  Cavity  Flows  Science  and  Engineering  Directorate 

Marshall  Space  Flight  Center 


- Inject  cooler  flow  near  inner-race  contact  areas  to  reduce  ball  temperatures 


20  April  1993 


- Three-dimensional  isothermal  analysis  using  FDNS 


- Effectiveness  determined  by  maximum  temperature  in  bearing  area 


Elastohvdrodynamic  film  (Q)  Lowered  shoulder  Coolant  injection 

None  (0.56  Btu/sec)  255  °R  21 3 °R 

Partial  (0.23  Btu/sec)  240  °R  204  °R 

Full  (0.028  Btu/sec)  231  °R  198  °R 
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OF  POOR  QUALITY 
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CFD  ANALYSIS  OF  ATD  BEARING  FLOWS 
FLOW  STREAMLINES  COLORED  BY  FLUID  TEMPERATURE 


N95-  23618 
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NAVIER-STOKES  FLOW  FIELD  ANALYSIS  OF  COMPRESSIBLE  FLOW  IN  A 
PRESSURE  RELIEF  VALVE 


f- 


i'i 


Bruce  Vu  & Ten-See  Wang 
Computational  Fluid  Dynamics  Branch 
NASA  - Marshall  Space  Flight  Center 


Ming-Hsin  Shih  & Bharat  Soni 

Engineering  Research  Center  for  Computational  Field  Simulation 
Mississippi  State  University 


Abstract 

The  present  study  was  motivated  to  analyze  the  complex  flow  field  involving  gaseous 
oxygen  (GOX)  flow  in  a relief  valve  (RV).  The  9391  RV,  pictured  in  Figure  1,  was  combined 
with  the  pilot  valve  to  regulate  the  actuation  pressure  of  the  main  valve  system.  During  a 
high-pressure  flow  test  at  Marshall  Space  Flight  Center  (MSFC)  the  valve  system  developed 
a resonance  chatter,  which  destroyed  most  of  the  valve  body.  Figures  2-4  show  the  valve  body 
before  and  after  accident.  It  was  understood  that  the  subject  RV  has  never  been  operated  at  5500 
psia.  In  order  to  fully  understand  the  flow  behavior  in  the  RV,  a computational  fluid  dynamics 
(CFD)  analysis  is  carried  out  to  investigate  the  side  load  across  the  piston  sleeve  and  the 
erosion  patterns  resulting  from  flow  distribution  around  piston/nozzle  interface. 


Grid  Topology 

The  safety  RV  consists  of  a main  cylinder  and  a piston,  with  a smaller  diameter  inlet. 
An  intersection  technique  was  developed  to  model  the  piston-cylinder  configuration  (Figure  5). 
To  simplify  the  geometry,  the  diameter  of  the  cylinder  is  kept  constant. 

An  O-type  grid  in  the  axial  plane  was  initially  considered  for  this  geometry. 
However,  it  would  become  very  difficult  to  generate  grid  lines  around  the  piston  on  the  upper 
part  of  the  main  cylinder,  if  not  impossible.  H-type  grid  was  then  chosen  to  model  this  internal 
flow  geometry.  The  main  cylinder  was  cut  into  half  at  the  plane  of  symmetry  to  reduce  the  size 
of  domain.  It  was  again  cut  into  halves  at  the  bottom  face  of  the  piston  to  divide  the 
computational  domain  into  upper  part  and  lower  part.  To  model  the  field  geometry,  the 
descritization  was  carried  out  into  five-block  zonal  grid;  the  inlet  itself  formed  a block,  the 
lower  part  of  the  main  cylinder  formed  another  block,  and  the  upper  part  of  the  main  cylinder 
was  cut  into  3 more  blocks.  The  5-block  grid  is  shown  in  Figures  6-7.  Compared  to  the  original 
O-type  grid,  this  H-type  grid  topology  greatly  reduced  the  grid  distortion,  especially  near  the 
piston  . 


Grid  Generation 

GENIE++  (Ref.  1-3),  a general  purpose  three-dimensional  grid  generation  package,  was 
used  to  generate  the  grid  for  thi6  geometry.  GENIE++  is  the  Mississippi  State  University 
updated  version  of  INGRID  (Ref.  4-5)  developed  by  Arnold  Engineering  Development  Center 
(AEDC). 

In  order  to  perform  the  surface  intersections  of  tire  piston,  as  welt  as  the  inlet,  with  the 
main  cylinder,  an  intersection  algorithm  with  Newton-Ralphson  method  was  used  to  obtain 
the  intersection  curves.  Weighted  transfinite  interpolation  (Weighted  TFl)  (Ref.  2)  algorithm 
is  used  to  generate  the  algebraic  grid.  Weighted  TFI  can  be  formulated  as  uniform  TFI  with  grid 
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distribution  mesh,  where  the  grid  distribution  mesh  is  obtained  by  performing  uniform  TFI  on 
normalized  arc  length  distribution  on  associated  boundaries  (or  surfaces  in  volume  grid). 

Since  the  selected  grid  topology  reduced  the  distortion  of  grid  lines,  the  resulting 
algebraic  grid  was  very  satisfactory,  and  no  elliptic  smoothing  was  performed  for  the  present 
computation.  However,  for  future  grid -dependent  study,  elliptic  solver  will  be  applied  to 
refine  the  local  grids  while  maintaining  a packed,  viscous  grid  on  the  surface. 

Governing  Equations  and  Computational  Scheme 

The  present  numerical  simulation  uses  a non-staggered  grid,  pressure  based  transport 
equation  solver  with  an  extended  version  of  two-equation  k-e  turbulence  model.  While  the 
computer  code  has  all-speed  capability  for  both  compressible  and  incompressible  flows,  the 
present  study  only  uses  the  compressible  feature.  The  basic  equations  employed  to  describe  the 
momentum  and  heat  transfer  in  the  computational  domain  are  the  three-dimensional  Reynolds- 
averaged  transport  equations.  To  solve  the  system  of  coupled  nonlinear  partial  differential 
equations,  it  uses  finite  difference  approximations  to  establish  a system  of  linearized  algebraic 
equations.  An  adaptive  up  winding  scheme  is  utilized  to  model  the  convective  terms  of  the 
momentum,  energy  and  continuity  equations,  which  is  based  on  the  second  and  fourth  order 
central  differencing  with  artificial  dissipation.  Discretization  of  viscous  fluxes  and  source 
terms  uses  a second-order  central  difference  approximation.  For  velocity-pressure  coupling,  the 
present  solution  procedure  employs  pressure-based,  predictor  followed  by  multi-corrector 
approach.  Details  of  the  present  numerical  methodology  are  given  by  Wang  and  Chen  (Ref.  6). 

Due  to  symmetry,  the  computational  domain  occupies  only  the  front  half  of  the  RV. 
Along  all  solid  walls,  no-slip  condition  is  applied  for  velocities,  and  temperature  is  assumed 
constant.  For  near-wall  turbulence  treatment,  it  uses  a wall  function  with  modified  flux  source 
and  a velocity  profile  capable  of  providing  a smooth  transition  between  logarithmic  law-of- 
wall  and  linear  variation  in  viscous  sublayer.  Such  a treatment  significantly  reduces  the  flux 
dependence  on  the  near-wall  spacing.  The  inlet  conditions  are  fully  developed  profiles  for 
velocities  and  turbulence  parameters,  and  the  outlet  conditions  satisfy  the  conservation  of  mass. 


Result  and  Discussion 

The  preliminary  computations  have  been  performed  to  simulate  the  flow  field  of  GOX 
in  the  9391  RV  at  5500  psia  and  1000°  R.  Results  indicated  no  viscous  heating  due  to  low 
temperature  gradients  near  the  piston  surface  (Figure  8).  The  surface  pressure  contours  in  Figure 
9 also  indicated  an  insignificant  side  load  across  the  piston  sleeve.  The  force  obtained  from 
integrating  all  pressure  points  around  the  piston  surface,  from  the  bottom  up  to  the  piston  sleeve 
is  found  to  be  only  70  lbf,  under  this  adverse  condition.  The  velocity  vectors,  magnitude,  and 
Mach  contours  are  shown  in  Figures  10-12,  respectively.  Finally,  the  vortex  formations  in 
Figure  13-14  predicts  reasonable  erosion  patterns  in  the  gap  between  the  cylinder  elbow  and  the 
bottom  of  the  piston.  Evidently  these  patterns  are  in  agreement  with  the  damaged  hardware 
which  indicates  clear  signs  of  burns  and  scratches  near  the  piston/ throat  region. 
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Figure  2.  Valve  photo  (before  accident) 
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Figure  3.  Side-view  photo  (after  accident) 
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Figure  5.  Intersection  technique 
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Figure  6.  Three-dimensional  grid  (side-view) 
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Figure  7.  Three-dimensional  grid  (end-view) 
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Figure  8.  Temperature  contours 
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Figure  9.  Pressure  contours 
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Figure  10.  Velocity  vectors 
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Figure  11.  Velocity  contours 
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Figure  12.  Mach  contours 
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Figure  13.  Particle  tracing 
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Figure  14.  Particle  tracing  (close-up) 
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COMPUTATIONAL  FLUID  DYNAMICS  ( CFD ) ANALYSES  IN  SUPPORT  OF  p_  ■ J f 
SPACE  SHUTTLE  MAIN  ENGINE  ( SSME ) HEAT  EXCHANGER  ( HX ) VANE 

CRACKING  INVESTIGATION 

R.  Garcia,  T.  Benjamin,  and  J.  Cornelison 
NASA/Marshall  Space  Flight  Center,  Alabama 

A.  J.  Fredmonski 

Pratt  & Whitney,  West  Palm  Beach,  Florida 


Integration  issues  involved  with  installing  the  alternate 
turbopump  (ATP)  High  Pressure  Oxygen  Turbopump  (HPOTP)  into  the 
SSME  have  raised  questions  regarding  the  flow  in  the  HPOTP 
turnaround  duct  (TAD).  Steady-state  Navier-Stokes  CFD  analyses 
have  been  performed  by  NASA  and  Pratt  & Whitney  (P&W)  to  address 
these  questions.  The  analyses  have  consisted  of  two-dimensional 
axisymmetric  calculations  done  at  Marshall  Space  Flight  Center 
and  three-dimensional  calculations  performed  at  P&W.  These 
analyses  have  identified  flowfield  differences  between  the 
baseline  ATP  and  the  Rocketdyne  configurations.  The  results 
show  that  the  baseline  ATP  configuration  represents  a more 
severe  environment  to  the  inner  HX  guide  vane.  This  vane  has 
limited  life  when  tested  in  conjunction  with  the  ATP  but 
infinite  life  when  tested  with  the  current  SSME  HPOTP.  The  CFD 
results  have  helped  interpret  test  results  and  have  been  used  to 
assess  proposed  redesigns.  This  paper  includes  details  of  the 
axisymmetric  model,  its  results,  and  its  contribution  towards 
resolving  the  problem. 
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Conclusion/future  work 
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- Differences  between  Rocketdyne  and  Pratt  & Whitney  configuration 

Identify  configurational  changes  to  the  turbopump  to  reduce  unsteady  loading 


SSME  POWERHEAD  COMPONENT  ARRANGEMENT 
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Hot  fire  testing  (Pratt  & Whitney  and  SSC) 

- Actual  environment,  duration  testing 

- Provides  HX  vane  strain  levels 


CD  W 


:s  ^ 

£ 5 

-a  u 

< I 


j=. 

fr  £ 

O TO 

0)  * E ® 

ifie 

«&“5  Ju 

■s»§§* 

«5  P .£ 

»«  E j 
.£=:£> 

"-1° 


2 ^ 
<9  XJ 

5 c 

<-  fT3 

U oi 


S &3D 

O »•  O T! 
re  2 E ‘3 

2 CJ  <J5  E 


cn: 

1-4, 

0/3 

><! 

hJ 

< 

Z 

< 

Q 

tu, 

U 
• • 
E 

U 

o! 

c4 

C4 

o< 

< 


<u 

a 

.35 

£ 


ns 

>-< 

Cn 


Q 

cn 

H 

U 

cd  < 
Oi  W 

H P* 


CD 

*4 

CD 

r— < 

rfl 

a 

03 

<u 

fi 
• 1^ 

x> 

u 

3 

4-» 

a» 

to 

rs 

44 

CD 

l 


3 

£ 

M 

a/ 

i— t 

3 

w 

”0 

0) 


ns 

O 


a» 

£ 

o 

a> 

60 

cu 

c 

• H 

'a; 

CD 

OS 

£> 


O 

£> 

O 


a> 


u 

> 

a> 


o 

3 

73 

73 

a> 

73 

•l-H 

> 

o 

u 

Oh 


<u 

3 

73 

«*-• 

0) 

p*s 

cu 

o 

g 

CD 

CD 

»— ( 
ns 
S3 
ns 

a> 

3 


B -e 

^ 5 


CD 

0/ 

34 

0 

4-> 

cn 

1 

a> 


3 .5 

03 

Z 

a> 
60 
ns 

-M 

CD 
0/ 


60 

3 

♦pH 

cn 


CD 

O 

• pH 
4-» 

cu 

• *-H 

73 

cu 

Vh 

a. 

Jh 

a> 

rH 

pi 

W 

73 

cu 


73 

O 

£ 

73 

cu 


Sh 

CU 

> 


cn 

a 

w 

-I 

Oh 

w 

P4 

u 

tin 

cn 


CD 

• *H 

CD 

p*H 

ns 

3 

f6 

Q 

h 

U 

73 

ns 

a» 

4-* 

CD 

<J 

• 4 
4-» 

a> 

a 

a 

£• 

•pH 

X 


CD 

S3 

O 

* n-4 
4— > 

ns 

j-( 

33 

60 

PP 

S3 

O 

cu 

cu 

c 

• I— I 

Id 

^ ss 

-m  ns 

JP  S3 

03  CD 

Oh  £ 

ns  k. 

cu  O 


73 

3 

3 

O 

J-H 

os 

3 


O 

4-* 

73 

0/ 

73 

* r-4 
> 
O 

Vh 


a/ 

JD 

CD 

0> 

U 

C 

cu 

S-H 

CU 


S3 

cu 

73 


73 

cu 

CD 


£ D 


CD 

S3 

O 

• *-H 
J— I 

03 

5 ^ 

bJD  5i 


S3 

o 

u 

6JD 


CD 

cu 

4-» 

73 

SP 

03 

V 

CD 

cu 

X 


73 

0) 

CD 

o 


73 


03 

cu 

-4-) 

03 

> 

cu 


73 

cu 

CD 

D 


rC 


ns 

j-i 

6 

CD 

CD 

ns 

S3 

ns 

P 

Mh 


ns 

45 

u 

CD 


2 

’aJ 

• 

M— < 

O 


S3 

O 

4-4 

- r,  2 

§■  U « 

2 ^ 72 

Oh  £P  2 


X 

03 


CU 

.a; 


£ cu 


a> 

£ 

t 

CD 

Q 

cn 


73 

CU 

73 

T> 

O 


< 

cn 

D 


a> 

S 

cu 

O 

g 

CD 

• *4 
CD 

2 

S3 

D 

Uh 

U 

73 

ns 

cu 

44 

CD 

fi 

3 

V 

u 

• hH 

M 

«*4 

<D 

a 

a 

CD 

• i^ 

X 

< 


106 


Used  to  identify  unsteady  flow  features  not  captured  by  the  steady  analysis 
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Code  solutions  used  to  obtain  relative  comparisons  between  configurations 

- Absolute  quantities  treated  with  caution 

Goal  was  to  try  to  match  the  Rocketdyne  velocity  and  turbulence  fields 
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- Velocity  and  turbulence  intensity  profile  at  the  turbopump-to-engine  interface 

- Velocity  and  turbulence  intensity  profile  at  the  HX  coils 
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To  assure  that  the  proposed  fix  is  not  impaired  functionally  by  3D  effects 
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channel  height  ( 1.  = o.d.  wall) 
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3D  FLOW  ANALYSIS  OF  THE  ALTERNATE  SSME  HPOT  TAD 


C.A.Kubinski 

Government  Engines  & Space  Propulsion  Division 
West  Palm  Beach,  Florida 


ABSTRACT 


This  paper  describes  the  results  of  numerical  flow  analyses  performed  in  support  of  design 
development  of  the  Space  Shuttle  Main  Engine  Alternate  High  Pressure  Oxidizer  Turbine 
Turn-around  duct  (TAD).  The  flow  domain  has  been  modeled  using  a 3D,  Navier -Stokes, 
general  purpose  flow  solver.  The  goal  of  this  effort  is  to  achieve  an  alternate  TAD  exit  flow 
distribution  which  closely  matches  that  of  the  baseline  configuration.  3D  Navier  Stokes  CFD 
analyses  were  employed  to  evaluate  numerous  candidate  geometry  modifications  to  the  TAD 
flowpath  in  order  to  achieve  this  goal.  The  design  iterations  are  summarized,  as  well  as  a de- 
scription of  the  computational  model,  numerical  results  and  the  conclusions  based  on  these 
calculations. 
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Workshop  for  Computational  Fluid  Dynamics  Applications  in 

Rocket  Propulsion,  NASA  MSFC 
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Cheryl  A.  Kubinski 
Pratt  & Whitney  (GESP) 
April  20,  1 993 
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SSME  ATD  HPOT  TAD  REDESIGN 
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SSME  ATP  HPOT  TAD  REDESIGN 

Problem  Overview 
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Unsteadiness/turbulence 

Turbine  thermal  profiles 


SSME  ATP  HPOT  TAD  REDESIGN 

TAD  Inlet  Gas  Profile  Differences 
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SSME  ATD  HPOT  TAD  REDESIGN 

Code  Verification 


Predicted  Static  Pressure  With  Trip  Strips 


SSME  ATD  HPOT  TAD  REDESIGN 

Code  Verification 


Measured  vs  Predicted  Nusselt  Number  For 
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REFLEQX  Axisymmetric  Analysis  (MSFC) 


rocketdyne  model. 


PROSTAR  2.1 


Baseline  ATD...converged 


SSME  ATP  HPOT  TAD  REDESIGN 

Differences  in  Velocity  Distributions  at  Interface  Plane 


REFLEQX  Axisymetric  Analysis  (MSFC) 


SSME  ATP  HPOT  TAD  REDESIGN 

Differences  in  Velocity  Distributions  at  Interface  Plane 
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3D  Analysis 


SSME  ATD  HPOT  TAD  REDESIGN 


SSME  ATP  HPOT  TAD  REDESIGN 

CFD  Models  Differ  in  Prediction  of  Flow  Separation  For  Modified  TAD 
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SSME  ATD  HPOT  TAD  REDESIGN 


PROSTAR  2.1 
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ATD  Redesign 


SSME  ATP  HPQT  TAD  REDESIGN 

Axisymmetric  Analysis  Indicates  Velocity  Near  Inner  Vane  Still  At  Goal  Level 
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Percent  Span  (ATD  Baseline) 


SSME  ATP  HPQT  TAD  REDESIGN 

3D  Analysis  ofATD  Redesign  TAD  Indicates  Velocity  Near  Inner  Vane  At  Goal 
Level 
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UNSTEADY  FLOW  SIMULATIONS  IN  SUPPORT  OF 
THE  SSME  HEX  TURNING  VANE  CRACKING 
INVESTIGATION  WITH  THE  ATD  HPOTP 

by 


N.  S.  Dougherty,  D.  W.  Burnette,  and  J.  B.  Holt 
Rockwell  International 
Space  Systems  Division 
Huntsville,  AL  35806 


and 


T.  Nesman 
MSFC/ED33 

Marshall  Space  Flight  Center,  AL  35812 


ABSTRACT 


Unsteady  flow  computations  are  being  performed  with  the  P&W  (ATD)  and  the 
Rocketdyne  baseline  configurations  of  the  SSME  LO2  turbine  turnaround  duct  (TAD)  and 
heat  exchanger  (HEX).  The  work  is  in  support  of  the  HEX  inner  turning  vane  cracking 
investigation.  Fatigue  cracking  has  occurred  during  hot  firings  with  the  P&W  configuration 
on  the  HEX  inner  vane,  and  it  appears  the  fix  will  involve  changes  to  the  TAD  splitter  vane 
position  and  to  the  TAD  inner  wall  curvature  to  reduce  the  dynamic  loading  on  the  inner 
vane.  Unsteady  flow  computations  on  the  P&W  baseline  and  fix  and  on  the  Rocketdyne 
baseline  reference  follow  steady-flow  screening  computations  done  by  MSFC/ED32  on 
several  trial  configurations  arriving  at  the  fix. 

The  P&W  TAD  inlet  velocity  profile  has  a strong  radial  velocity  component  that  directs  the 
flow  toward  the  inner  wall  and  raises  the  local  velocity  a factor  of  two  and  the  dynamic 
pressure  a factor ,of  four.  The  fix  is  intended  to  redistribute  the  flow  more  evenly  across 
the  HEX  inner  and  outer  vanes  like  the  Rocketdyne  baseline  reference.  Vane  buffeting  at 
frequencies  around  4,000  Hz  is  the  leading  suspected  cause  of  the  problem.  Our 
simulations  (work  in  progress)  are  being  done  with  the  USA  2D  axisymmetric  code 
approximating  the  flow  as  axisymmetric  u+v  2D  (axial,  u,  and  radial,  v,  components  only). 
The  HEX  coils  are  included  in  the  model  to  make  sure  the  fix  does  not  adversely  affect  the 
HEX  environment. 

Turbulent  kinetic  energy,  k,  levels  where  k = 1/2  v’  rms2  are  locally  as  high  as  10,000 
ft2/sec2  for  the  P&W  baseline  at  the  engine  interface  (between  the  TAD  and  HEX)  at  the 
HEX  inner  vane  location.  However,  k is  less  than  8,000  on  the  HEX  outer  vane  and  only 
about  4,500  on  the  HEX  inner  vane  for  the  Rocketdyne  baseline.  Unsteady  turbulence 
intensity,  v’rms/v,  and  pressure,  p’,  are  being  computed  in  the  present  computations  to 
compare  with  steady-flow  Reynolds-averaged  computations  where  p’rms  = const  (pk)  for 
overall  rms  random  turbulence  from  0.1  to  12,000  Hz  frequency.  Random  overall  static, 
p'rms  fluctuations  as  large  as  1.7  psi  are  estimated  from  k on  the  HEX  inner  vane  for  the 
P&W  baseline  configuration  but  only  about  0.7  psi  for  the  Rocketdyne  configuration. 
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Rockwell  ENGINE  HEX  TURNING  VANE  CRACKING  INVESTIGATION 
International UNSTEADY  CFD  ANALYSIS 

Space  Systems  Division 
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TURBULENCE  TREATMENT  IS  ALWAYS  A QUESTION 
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- LDV  Measurements  Obtained  From  Rocketdyne  and  NASA/MSFC 
at  Q/Qd=0.8,  1.0,  and  1.15 

- Non-Dimensional  Velocities  (Vel/Utip) 
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Three-Dimensional  Flow  Analysis 
Inside  Consortium  Impeller 
at  Design  and  Off-Design  Conditions 


C.Hah,  J.  Loellbach,  F-L.  Tsung,  and  D.  A.  Greenwald 
NASA  Lewis  Research  Center 

R.  Garcia 

NASA  Marshall  Space  Flight  Center 


Three-dimensional  flow  fields  inside  the  Consortium  impeller 
were  analyzed  with  a Navier-Stokes  code.  The  numerical  results 
at  the  design  and  off-design  conditions  are  compared  with  the 
experimental  data. 


195 


W 

Q 

>—< 

2 Q 
& £ 
1/5  ^ 


Z 

O 

i/5 

Q 


c/3 

<! 

Z 

< 

S tf  2 

j H c 

^ »-J  |-h 

5 a £ 

O 5 Z 
55  B O 

Z s U 

■Ed  Z _ 

S P Z 

W 0 

2H  w 

Q js5  w 

wog 

HM® 

« Z EK 
aofe 
H U O 


Z 

Ed 

Ed 

fig 

O 


Q 


Q 

Z 

< 

«r» 

a 

z 

3 

co 

r 


< 

S3 

J 

J 

Ed 

O 

J 

• 

*-5 


< 

X 

u 


OS 

cd 

H 

z 

cd 

u 


u 


< 

Ed 

C/3 

Ed 


5/3 

£ 

J 

< 

C/3 

< 

Z 


< 

M 

U 

Csi 

< 


X 

Ed 

H 

Z 

Ed 

U 

H 

X 

a 

E=1 

Ed 

U 

< 

A 

c/3 

J 

J 

< 

co 

< 


c 

C/3 

< 

z 


196 


Objective 


198 


Computational  grid 


■© 

c 

C3 


C/3 

Q 

a 


© 


fS 

X 

rr 

<N 

© 


<N 


X 

Tf 

<N 

© 

© 


199 


Computational  Grid 


Flow  Cases  Studied 
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Nondimensional  Radial  Velocity 
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OPTIMIZATION  OF  A CENTRIFUGAL  IMPELLER  DESIGN 

THROUGH  CFD  ANALYSIS  ^ * 

W.  C.  Chen,  A.  H.  Eastland,  D.  C.  Chan  , / 

Rockwell  International,  Rocketdyne  Division  V 

R.  Garcia  ■) 

NASA,  Marshall  Space  Flight  Center  ; ' 

This  paper  discusses  the  procedure,  approach  and  Rocketdyne  CFD  results  for 
the  optimization  of  the  NASA  consortium  impeller  design.  Two  different 
approaches  have  been  investigated.  The  first  one  is  to  use  a tandem  blade 
arrangement,  the  main  impeller  blade  is  split  into  two  separate  rows  with  the 
second  blade  row  offset  circumferentially  with  respect  to  the  first  row.  The 
second  approach  is  to  control  the  high  losses  related  to  secondary  flows  within 
the  impeller  passage.  Many  key  parameters  have  been  identified  and  each 
consortium  team  member  involved  will  optimize  a specific  parameter  using  3-D 
CFD  analysis.  Rocketdyne  has  provided  a series  of  CFD  grids  for  the  consortium 
team  members.  SECA  will  complete  the  tandem  blade  study,  SRA  will  study  the 
effect  of  the  splitter  blade  solidity  change,  NASA  LeRC  will  evaluate  the  effect  of 
circumferential  position  of  the  splitter  blade,  VPI  will  work  on  the  hub  to  shroud 
blade  loading  distribution,  NASA  Ames  will  examine  the  impeller  discharge 
leakage  flow  impacts  and  Rocketdyne  will  continue  to  work  on  the  meridional 
contour  and  the  blade  leading  to  trailing  edge  work  distribution.  This  paper  will 
also  present  Rocketdyne  results  from  the  tandem  blade  study  and  from  the 
blade  loading  distribution  study.  It  is  the  ultimate  goal  of  this  consortium  team  to 
integrate  the  available  CFD  analysis  to  design  an  advanced  technology  impeller 
that  is  suitable  for  use  in  the  NASA  Space  Transportation  Main  Engine  (STME) 
fuel  turbopump. 
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USE  OF  BLADE  LEAN  IN  TURBOMACHINERY  REDESIGN 


John  Moore,  Joan  G.  Moore,  and  Alex  Lupi 


Mechanical  Engineering  Department 
Virginia  Polytechnic  Institute  and  State  University 
Blacksburg,  Virginia  24061-0238 


Blade  lean  is  used  to  improve  the  uniformity  of  exit  flow  distributions  from 
turbomachinery  blading.  In  turbines,  it  has  been  used  to  control  secondary  flows  by  tailoring 
blade  turning  to  reduce  flow  overturning  and  underturning  and  to  create  more  uniform  loss 
distributions  from  hub  to  shroud. 

In  the  present  study,  the  Pump  Consortium  centrifugal  impeller  has  been  redesigned  using 
blade  lean.  The  flow  at  the  exit  of  the  baseline  impeller  had  large  blade-to-blade  variations, 
creating  a highly  unsteady  flow  for  the  downstream  diffuser.  Blade  lean  is  used  to  redesign  the 
flow  to  move  the  high  loss  fluid  from  the  suction  side  to  the  hub,  significantly  reducing  blade-to- 
blade  variations  at  the  exit. 


Axial  Flow  Turbine  Stators 
Consortium  Pump  Impeller  Problem 
Secondary  Flow  Analysis  for  a Rotor 
Stable  Location  of  High  Loss  Fluid 
Impeller  Redesign 
Improved  Performance 


Use  of  Blade  Lean 
in 

Axial  Flow  Turbine  Stators 
VP  in  Cross-Sections 
Controlling  Exit  Loss  Distributions 
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Expect  more  losses  Opposing  radially  Spreading  loss  spanwise 

in  hub/suction  side  inward  flow  in 

corner  region.  suction  side  boundary  layer. 
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MODEL  MB  LOW  SOLIDITY  VANES 
BELLMOUTH  END  WALL  (Cl) 
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Effect  of  compound  lean  on  stator 


PRESSURE  A A SUCTION 
SURFACE  V V SURFACE 


Phase  1/ 

Model  D stator 


CONSORTIUM  IMPELLER 
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Equations  for  Incompressible  Flow 
in  a Rotating  System 


Reduced  static  pressure 

pr  = p - I pa2 r2. 


Rotary  stagnation  pressure 

P*  = P + ^pW2  - Ipw2 r2 

Absolute  vorticity 

Q = VxV  = VxW  + 2uj 

Momentum,  inviscid  flow 

(\A/-V)W  + 2inx)H  = --Vpr 

P 
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Determining  the  Stable  Orientation  Vector 
for  Secondary  Vorticity  Suppression 
in  Rotating  Systems 

Secondary  Circulation,  Hawthorne 


a 

Os' 

2 

r i ap* 

+ w ap*' 

as 

w 

b.  > 

pW2 

Rn  3b 

< 

w dz 

From  momentum 

w x 0 = 1 Vp* 
P 


M ■ Vp*  = 0,  Vp*  1 W 


Generation  of  secondary  circulation  = 0 when 


W • 


l--VPr 

p 

V 


wxW]  x Vp* 


l.e.  the  component  of  the  vector 

Vpr  - ^ x W 
P 


perpendicular  to  the  relative  velocity 
points  to  the  stable  location  of  high  loss  fluid. 
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Consortium  Impeller 
Contours  of  P*  at  the  impeller  exit 
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DIFFUSER  VANE  EXCITATION  PARAMETER 
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Impeller  exit  Diffuser  inlet 
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Reduced  diffuser  vane  excitation  forces 
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CFD  PARAMETRIC  STUDY  OF  CONSORTIUM  IMPELLER 

Gary  C.  Cheng*,  Y.S.  Chen*,  R.  Garcia*,  and  R.W.  Williams* 


Abstract 


Current  design  of  high  performance  turbopumps  for  rocket  engines  requires  effective  and 
robust  analytical  tools  to  provide  design  impact  in  a productive  manner.  The  main  goal  of  this 
study  is  to  develop  a robust  and  effective  computational  fluid  dynamics  (CFD)  pump  model  for 
general  turbopump  design  and  analysis  applications.  A Finite  Difference  Navier-£tokes  flow 
solver,  FDNS,  which  includes  the  extended  k-£  turbulence  model  and  appropriate  moving 
interface  boundary  conditions,  was  developed  to  analyze  turbulent  flows  in  turbomachinery 
devices.  A second-order  central  difference  scheme  plus  adaptive  dissipation  terms  was  employed 
in  the  FDNS  code,  along  with  a predictor  plus  multi-corrector  pressure-based  solution  procedure. 
The  multi-zone,  multi-block  capability  allows  the  FDNS  code  to  efficiently  solve  flow  fields  with 
complicated  geometry.  The  FDNS  code  has  been  benchmarked  by  analyzing  the  pump 
consortium  inducer,  and  it  provided  satisfactory  results.  In  the  present  study,  a CFD  parametric 
study  of  the  pump  consortium  impeller  was  conducted  using  the  FDNS  code.  The  pump 
consortium  impeller,  with  partial  blades,  is  a new  design  concept  of  the  advanced  rocket  engines. 
The  parametric  study  was  to  analyze  the  baseline  design  of  the  consortium  impeller  and  its 
modification  which  utilizes  TANDEM  blades.  In  the  present  study,  the  TANDEM  blade 
configuration  of  the  consortium  impeller  considers  cut  full  blades  for  about  one  quarter  chord 
length  from  the  leading  edge  and  clocks  the  leading  edge  portion  with  an  angle  of  7.5  or  22.5 
degrees.  The  purpose  of  the  present  study  is  to  investigate  the  effect  and  trend  of  the  TANDEM 
blade  modification  and  provide  the  result  as  a design  guideline.  A 3-D  flow  analysis,  with  a 103 
x 23  x 30  mesh  grid  system  and  with  the  inlet  flow  conditions  measured  by  Rocketdyne,  was 
performed  for  the  baseline  consortium  impeller.  The  numerical  result  shows  that  the  mass  flow 
rate  splits  through  various  blade  passages  are  relatively  uniform.  Due  to  the  complexity  of  blade 
geometries,  the  TANDEM  blade  configurations  were  analyzed  with  the  multi-zone  grid  structure. 
Both  the  7.5°-  and  the  22. 5 “-clocking  TANDEM  blade  cases  utilized  a 80K  mesh  system.  The 
numerical  result  of  two  TANDEM  blade  modifications  indicates  the  efficiency  and  the  head  are 
worse  than  those  of  the  baseline  case  due  to  larger  flow  distortion.  The  gap  between  the 
TANDEM  blade  and  the  full  blade  allows  the  flow  passes  through  and  heavily  loads  the  pressure 
side  of  the  partial  blade  such  that  flow  reversal  occurs  near  the  suction  side  of  the  splitter.  The 
flow  split  at  the  exit  of  impeller  blades  is  very  non-uniform  for  TANDEM  blade  cases,  and  this 
will  greatly  induce  the  side  load  on  the  diffuser.  Therefore,  the  TANDEM  blade  modification 
in  the  present  CFD  analysis  does  not  improve  the  performance  of  the  consortium  impeller. 


SEC  A,  Inc.,  3313  Bob  Wallace  Ave.,  Suite  202,  Huntsville,  AL 
Engineering  Sciences,  Inc.,  4920  Corporate  Dr.,  Suite  K,  Huntsville,  AL 
ED  32,  NASA/Marshall  Space  Flight  Center,  Huntsville,  AL 
ED  32,  NASA/Marshall  Space  Flight  Center,  Huntsville,  AL 
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ELEVENTH  WORKSHOP  FOR  CFD  APPLICATIONS  IN  ROCKET  PROPULSION 

NASA/MSFC,  APRIL  20-22,  1993 


OBJECTIVE 
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TEST  CONFIGURATION  SETUP 
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STREAMWISE  DIRECTION 


BLADE-TO-BLADE  DIRECTION  (PRESSURE  TO  SUCTION) 
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NUMERICAL  METHOD 
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STANDARD  AND  EXTENDED  K -e  TURBULENCE  MODELS 
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The  Mesh  System  Layout  for  the  Consortium  Impeller  (Hub-to-Tip) 
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CYCLIC  TANDEM 

B.C.  BLADE  CYCLIC  B.C.  CYCLIC  B.C. 
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The  Mesh  System  Layout  for  22.5°  Clocking  TANDEM  Blade 
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Suction  Side  of  Full  Blade  Pressure  Side  of  Partial  Blade 
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Suction  Side  of  Partial  Blade  Pressure  Side  of  Full  Blade 


BASELINE 
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Velocity  Vectors  Near  the  Hub 
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Pressure  Side  of  Partial 


7.5°  TANDEM  Blade 


Velocity  Vectors  at  the  Mid  Span 


7.5°  TANDEM  Blade 
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Velocity  Vectors  Near  the  Shroud 


22.5°  TANDEM  Blade 


Velocity  Vectors  Near  the  Hub 


22.5°  TANDEM  Blade 


Velocity  Vectors  Near  the  Shroud 
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Velocity  Vectors  at  the  Exit  of  Impeller  Blade 
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N95- 23627 


Abstract  of  a proposed  presentation  at  workshop  for 
CFD  Applications  in  Rocket  Propulsion  to  be  held  at 
NASA  Marshall  Space  Flight  Center,  AL,  April  20-22,  1993 

INCOMPRESSIBLE  NAVIER-STOKES  CALCULATIONS  IN  PUMP  FLOWS 

Cetin  Kiris,  Leon  Chang 
MCAT  Institute,  Moffett  Field,  CA 

and 

Dochan  Kwak 

NASA-Ames  Research  Center,  Moffett  Field,  CA 

Flow  through  pump  components,  such  as  the  SSME-HPFTP  Impeller  and  an 
advanced  rocket  pump  impeller,  is  efficiently  simulated  by  solving  the  incompressible 
Navier-Stokes  equations.  The  solution  method  is  based  on  the  pseudocompressibility 
approach  and  uses  an  implicit-upwind  differencing  scheme  together  with  the  Gauss- 
Seidel  line  relaxation  method.  The  equations  are  solved  in  steadily  rotating  reference 
frames  and  the  centrifugal  force  and  the  Coriolis  force  are  added  to  the  equation  of  mo- 
tion. Current  computations  use  one-equation  Baldwin-Barth  turbulence  model  which 
is  derived  from  a simplified  form  of  the  standard  k - e model  equations.  The  resulting 
computer  code  is  applied  to  the  flow  analysis  inside  an  11-inch  SSME  High  Pressure 
Fuel  Turbopump  impeller,  and  an  advanced  rocket  pump  impeller.  Numerical  results 
of  SSME-HPFTP  impeller  flow  are  compared  with  experimental  measuremnts.  In  the 

I 

advanced  pump  impeller,  the  effects  of  exit  and  shroud  cavities  are  investigated.  Flow 
analyses  at  design  conditions  will  be  presented. 
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Summary 


Introduction 
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> Coarse/ medium  size  grid  for  engineering  purposes  (150K  - 600K). 
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converges  to  a steady  state 
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SSME-HPFTP  Impeller  Computations 
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,SME  HPFTP  Impeller  Grid 
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SSME-HPFTP  Impeller  Computations 
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Advanced  Rocket  Pump  ImpeUer 

Shrouded  impeller 

6 full  blades 
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Case  3 - Exit  cavities  + shroud  cavity 


Advanced  Impeller  Computations 
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EFFICIENCY  vs  RELATIVE  X 
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RELATIVE  X FROM  SHROUD 


MERIDIONAL  VELOCITY  @X=0.5  vs  RELATIVE  X 
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RELATIVE  ANGLE  FROM  SUCTION  SIDE 


Advanced  Impeller  Concept 
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Advanced  Impeller  Concept 
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no  cavity  (slip  b.c.) 


Advanced  Impeller  Concept 


^5 


a?"  QQQOQ3332222nn- 
t \ QQQQaQOQQOSSSnn 

Q W SV-r--cNjrvjc^(Nincnnn'T'T’T 


_3oooooggoa 
-^aoocsooaoa 
- ji  am  □ m o 1/1  o in  □ 

-MP^ONlflNONin 
' lo  in  ujujioifir«-p^n- 


aoaaa^^aooaaao- 


— -.  ~ n a 


oooaaaoa 


326 


OniG^'iAL  PAGE  IS 
OF  POOR  QUALITY 


MERIDIONAL  VELOCITY  @X=0.3  vs  RELATIVE  ANGLE 
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RELATIVE  ANGLE  FROM  SUCTION  SIDE 


MERIDIONAL  VELOCITY  <5)X=0.5  vs  RELATIVE  ANGLE 
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RELATIVE  ANGLE  FROM  SUCTION  SIDE 


MERIDIONAL  VELOCITY  @X=0.7  vs  RELATIVE  ANGLE 
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RELATIVE  ANGLE  FROM  SUCTION  SIDE 


EFFICIENCY  vs  RELATIVE  X 
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RELATIVE  X FROM  SHROUD 


MERIDIONAL  VELOCITY  vs  RELATIVE  X 
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RELATIVE  X FROM  SHROUD 


RELATIVE  FLOW  ANGLE  vs  X 
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RELATIVE  X FROM  SHROUD 


Advanced  Impeller  Concept 
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Advanced  Impeller  Concept 
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no  cavity  ( slip  b.c.)  with  exit  cavities 


Advanced  Impeller  Concept 
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with  exit  cavities 


Advanced  Impeller  Concept 
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Advanced  Impeller  Concept 


337 


ORIGINAL  PAGE  IS 
OF  POOR  QUALITY 


with  exit  cavities 


Solution  procedure  for  rocket  engine  pump  analysis  was  validated 
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THE  INFLUENCE  OF  SWIRL  BRAKES  AND  A TIP  DISCHARGE  ORIFICE  ON 
THE  ROTORDYNAMIC  FORCES  GENERATED  BY  DISCHARGE-TO-SUCTION 
LEAKAGE  FLOWS  IN  SHROUDED  CENTRIFUGAL  PUMPS 


J.  M.  Sivo,  A.  J.  Acosta,  C.  E.  Brennen,  T.  K.  Caughey 

California  Institute  of  Technology 
Division  of  Engineering  and  Applied  Science 
Pasadena,  California  91125 

ABSTRACT 

Recent  experiments  conducted  in  the  Rotor  Force  Test  Facility  at  the  California  Insti- 
tute of  Technology  have  examined  the  effects  of  a tip  leakage  restriction  and  swirl  brakes 
on  the  rotordynamic  forces  due  to  leakage  flows  on  an  impeller  undergoing  a prescribed 
circular  whirl.  The  experiments  simulate  the  leakage  flow  conditions  and  geometry  of  the 
Alternate  Turbopump  Design  (ATD)  of  the  Space  Shuttle  High  Pressure  Oxygen  Turbop- 
ump and  are  critical  to  evaluating  the  pump’s  rotordynamic  instability  problems. 

Previous  experimental  and  analytical  results  have  shown  that  discharge-to-suction  leak- 
age flows  in  the  annulus  of  a shrouded  centrifugal  pump  contribute  substantially  to  the 
fluid  induced  rotordynamic  forces.  Also,  previous  experiments  have  shown  that  leakage 
inlet  (pump  discharge)  swirl  can  increase  the  cross-coupled  stiffness  coefficient  and  hence 
increase  the  range  of  positive  whirl  for  which  the  tangential  force  is  destabilizing.  In  recent 
experimental  work,  the  present  authors  demonstrated  that  when  the  swirl  velotity  within 
the  leakage  path  is  reduced  by  the  introduction  of  ribs  or  swirl  brakes,  then  a substantial 
decrease  in  both  the  destabilizing  normal  and  tangential  forces  could  be  achieved. 

Motivation  for  the  present  research  is  that  previous  experiments  have  shown  that  re- 
strictions such  as  wear  rings  or  orifices  at  pump  inlets  affect  the  leakage  forces.  Recent 
pump  designs  such  as  the  Space  Shuttle  Alternate  Turbopump  Design  (ATD)  utilize  tip 
orifices  at  discharge  for  the  purpose  of  establishing  axial  thrust  balance.  The  ATD  has 
experienced  rotordynamic  instability  problems  and  one  may  surmise  that  these  tip  dis- 
charge orifices  may  also  have  an  important  effect  on  the  normal  and  tangential  forces  in 
the  plane  of  impeller  rotation.  The  present  study  determines  if  such  tip  leakage  restrictions 
contribute  to  undesirable  rotordynamic  forces. 

Additional  motivation  for  the  present  study  is  that  the  widening  of  the  leakage  path 
annular  clearance  and  the  installation  of  swirl  brakes  in  the  ATD  has  been  proposed  to  solve 
its  instability  problems.  The  present  study  assesses  the  effect  of  such  a design  modification 
on  the  rotordynamic  forces. 

The  experimental  apparatus  consists  of  a solid  or  dummy  impeller,  a housing  instru- 
mented for  pressure  measurements,  a rotating  dynamometer  and  an  eccentric  whirl  mech- 
anism. The  solid  impeller  is  used  so  that  leakage  flow  contributions  to  the  forces  are 
measured,  but  the  main  throughflow  contributions  are  not  experienced.  The  inner  surface 
of  the  housing  has  been  modified  to  accommodate  meridional  ribs  or  swirl  brakes  within 
the  leakage  annulus.  In  addition,  the  housing  has  been  modified  to  accommodate  a dis- 
charge orifice  that  qualitatively  simulates  one  side  of  the  balance  piston  orifice  of  the  Space 
Shuttle  ATD. 

Results  indicate  the  detrimental  effects  of  a discharge  orifice  and  the  beneficial  effects 
of  brakes.  Plots  of  the  tangential  and  normal  forces  versus  whirl  ratio  show  a substantial 
increase  in  these  forces  along  with  destabilizing  resonances  at  some  positive  whirl  ratios 
when  a discharge  orifice  is  added.  When  brakes  are  added,  some  of  the  detrimental  effects  of 
the  orifice  are  reduced.  For  the  tangential  force,  a plot  versus  whirl  ratio  shows  a significant 
reduction  and  a destabilizing  resonance  appears  to  be  eliminated.  For  the  normal  force, 
although  the  overall  force  is  not  reduced,  again  a destabilizing  resonance  appears  to  be 
eliminated. 
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Single-auction  impeller  with  a balancing  chamber  on  the  back. 
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ATD  - HIGH  PRESSURE  OXIDIZER  TURBOPUMP 
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ROTORDYNAMIC  FORCES 


For  a circular  whirl  orbit: 

x*(t)  = ecos(Qt ) 
y*(t ) = esin(Qt ) 

F, \(t)  = \(Alx  + A’y!,)e 

Fl'(t)  = ±(-A'xy  + A;i)e 
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ROTORDYNAMIC  COEFFICIENTS 


M = Direct  Added  Mass 
C = Direct  Damping 
c = Cross-coupled  Damping 
K = Direct  Stiffness 
k = Cross-coupled  Stiffness 

k/C  = Whirl  Ratio 
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RING  - 
VOLUTE 
ADAPTER 


SEAL  RING 


ROTATING  SHROUD 


CALTECH  LEAKAGE  FORCE  TEST  APPARATUS 


TEST  MATRIX 


Table  1.  Tests  Without  Inlet  Swirl 


RPM  Q/u  Brahes  Q (GPM)  4> 


1000  -0.9  to  +0.9  0 0 0.0 

10  0.026 

20  0.052 

30  0.077 

4 0 0.0 

10  0.026 

20  0.052 

30  0.077 

8 0 0.0 

10  0.026 

20  0.052 

30  0.077 

2000  -0.6  to  +0.7  0 0 0.0 

10  0.013 

20  0.026 

30  0.039 

4 0 0.0 

10  0.013 

20  0.026 

30  0.039 

8 0 0.0 

10  0.013 

20  0.026 

30  0.039 
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B OGPM 
♦ 10GPM 
a 30  GPM 


QJ  0) 


a. 


a OGPM 
• 10  GPM 

■ 30  GPM 


ft/CO 


Figure  ( I ) Dimensionless  normal  and  tangential  forces  at  2000  RPM  with 
0 swirl  brakes  and  flow  rates  of  0,  10  and  30  GPM. 
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0/a> 


Figure  ( 2 ) Dimensionless  normal  and  tangential  forces  at  2000  RPM  with 
4 swirl  brakes  and  flow  rates  of  0,  10  and  30  GPM. 
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a o Brakes 
♦ 4 Brakes 

■ 8 Brakes 


b o Brakes 
• 4 Brakes 
■ 8 Brakes 


Figure  ( 3 ) Dimensionless  normal  and  tangential  forces  at  2000  RPM  and  a 
flow  rate  of  10  GPM  for  0,  4 and  8 swirl  brakes. 
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Figure  (4)  Dimensionless  normal  and  tangential  forces  at  2000  ft  PM  and  a 
flow  rate  of  30  GPM  for  0,  4 and  8 swirl  brakes. 
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CONCLUSIONS  FOR  PHASE  1 


1.  The  addition  of  brakes  reduces  the  destabilizing 
normal  force  for  all  flow  rates  tested. 

2.  For  flow  rates  below  $ = 0.025,  the  addition  of  brakes 
reduces  the  tangential  force  and  whirl  ratio. 
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ATD  - HIGH  PRESSURE  OXIDIZER  TURBOPUMP 
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RADIAL  TIP 
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CONCLUSIONS  FOR  PHASE  2 


1.  A tip  discharge  orifice  of  the  type  used  for  the 
Alternate  Turbopump  Design  (ATD)  of  the  Space 
Shuttle  High  Pressure  Oxygen  Turbopump  is  desta- 
bilizing. 

2.  The  design  modification  of  widening  the  leakage 
path  annular  clearance  and  installation  of  11  swirl 
brakes  in  the  ATD  would  reduce  some  of  the  detri- 
mental effects  of  the  orifice. 
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ADAPTATION  OF  THE  ADVANCED  SPRAY  COMBUSTION  CODE  f]  // 
TO  CAVITATING  FLOW  PROBLEMS  r- 

Pak-Yan  Liang 

Rocketdyne  Division,  Rockwell  International 


ABSTRACT 

A very  important  consideration  in  turbopump  design  is  the  prediction  and 
prevention  of  cavitation.  Thus  far  conventional  CFD  codes  have  not  been 
generally  applicable  to  the  treatment  of  cavitating  flows.  Taking  advantage  of  its 
two-phase  capability,  the  Advanced  Spray  Combustion  Code  is  being  modified 
to  handle  flows  with  transient  as  well  as  steady-state  cavitation  bubbles.  The 
volume-of-fluid  approach  incorporated  into  the  code  is  extended  and 
augmented  with  a liquid  phase  energy  equation  and  a simple  evaporation 
model.  The  strategy  adopted  also  successfully  deals  with  the  cavity  closure 
issue.  Simple  test  cases  will  be  presented  and  remaining  technical  challenges 
will  be  discussed. 
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VOLUME-OF-FLUID  TWO-PHASE  TRACKING 
SCHEME  IMPLEMENTED  IN  BOTH 
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SUMMARY  OF  GOVERNING  EQUATIONS 
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OVERALL  FLOW  CHART  FOR  ASCOMB 
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KEY  SUBROUTINES 
MODINP 
CALCUV 

CALCP 

OUTBC 

CALCF 


CALCENL 


CALCSC  (ITE) 
CALCSC  (IED) 

MODV1S 


CHEM 

CALCSP 

OUTBC 

CALCSC  (IEN) 
TEMPER 


MODPRO 


PRINT 
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ADVANTAGES  OF  VOF-APPROACH  OVER  OTHER 

CAVITATION  MODELS 
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• HIGHER  ORDER  EFFECTS  (DIFFERENT  p» /p g RATIOS,  FINITE 
EVAPORATION  RATES,  RECONDENSATION,  ROUGHNESS  EFFECTS 
ON  CAVITATION  INCEPTION  ETC.)  CAN  BE  ACCOMMODATED 


TYPICAL  CONVERGENCE  CHARACTERISTICS 
FOR  NONREACTING  CAVITATING  FLOW 
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STATUS  OF  THE  CAVITATION  UPGRADE  SUBTASK 
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LIQUID  EVAPORATION  TERM 


CONCLUDING  REMARKS 


Cavitation  Modeling  in  Euler  and  Navier-Stokes  Codes 


Manish  Deshpande,  Jinzhang  Feng,  and  Charles  L.  Merkle 
Propulsion  Engineering  Research  Center 
The  Pennsylvania  State  University 
State  College,  PA 


Many  previous  researchers  have  modeled  sheet  cavitation  by  means  of  a 
constant  pressure  solution  in  the  cavity  region  coupled  with  a velocity  potential 
formulation  for  the  outer  flow.  The  present  paper  discusses  the  issues  involved  in 
extending  these  cavitation  models  to  Euler  or  Navier-Stokes  codes.  The  approach 
taken  is  to  start  from  a velocity  potential  model  to  ensure  our  results  are 
compatible  with  those  of  previous  researchers  and  available  experimental  data, 
and  then  to  implement  this  model  in  both  Euler  and  Navier-Stokes  codes.  The 
model  is  then  augmented  in  the  Navier-Stokes  code  by  the  inclusion  of  the  energy 
equation  which  allows  the  effect  of  subcooling  in  the  vicinity  of  the  cavity  interface 
to  be  modeled  to  take  into  account  the  experimentally  observed  reduction  in  cavity 
pressures  that  occurs  in  cryogenic  fluids  such  as  liquid  hydrogen.  Although  our 
goal  is  to  assess  the  practicality  of  implementing  these  cavitation  models  in 
existing  three-dimensional,  turbomachinery  codes,  the  emphasis  in  the  present 
paper  will  center  on  two-dimensional  computations,  most  specifically  isolated 
airfoils  and  cascades.  Comparisons  between  velocity  potential,  Euler  and  Navier- 
Stokes  implementations  indicate  they  all  produce  consistent  predictions. 
Comparisons  with  experimental  results  also  indicate  that  the  predictions  are 
qualitatively  correct  and  give  a reasonable  first  estimate  of  sheet  cavitation  effects 
in  both  cryogenic  and  non-cryogenic  fluids.  The  impact  on  CPU  time  and  the  code 
modifications  required  suggests  that  these  models  are  appropriate  for 
incorporation  in  current  generation  turbomachinery  codes. 
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- Navier-Stokes  + Energy 
- Thermal  Effects 


Governing  Equations 
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Numerical  Scheme 


splitting. 


Euler  Analysis 
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- Use  Closure  Condition  From  Potential  Methods 


384 


</) 


C 0) 

o o 


a 

(/> 


a> 

> 

O 

385 


£ 

o 


o 

CD 


i 


NonLinear"  - B.C.'s  applied  on  Interface 

Computational  Domain  evolves 
with  solution 
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Airfoil 


Code  Validation 


Stokes/Euler  Comparison 


Pressure  Distribution 
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Cavity  Length  Comparison 
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Midchord  Cavitation 
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Effect  of  Cascade  Spacing 
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Cavity  Profile 


Thermodynamic  Boundary  Condition 
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Temperature  Depression  - Nitrogen 
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Cavity  Length  (L/c) 


Impact  of  Thermal  Effects  on  Cavitation 


Fraction  Chord  x/C 


Summary 
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Predictions  of  Panel,  Euler  and  N-S  are  similar 

- BL's  and  Turbulence  have  little  effect  on  Model 
Predictions 
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Model  predicts  magnitude  of  temperature  depression 
satisfactorily 


N95-  23631 


WORKSHOP  FOR  CFD  APPLICATIONS  IN  ROCKET 
PROPULSION  AT  NASA-MSFC 

Mehtab  M.  Pervaiz 
April  20,  1993 


Title:  An  Inducer  CFD  Solution  and  Effects  Associated  With  Cavitation 

This  presentation  describes  a CFD  analysis  for  an  Alternate  Turbopump  Devel- 
opment (ATD)  configuration.  The  analysis  consists  of  a coupled  configuration  of  the 
inducer  and  impeller.  The  work  presented  here  is  a joint  collaboration  of  J.  Garrett,  J. 
Kuryla  and  myself. 


Outline: 

This  view  graph  provides  an  outline  of  the  current  presentation.  I will  start  with  the 
ATD  configuration  for  the  inducer  and  impeller  and  the  corresponding  CFD  solution. 
Subsequently,  I will  describe  the  current  cavitation  modeling  approach  that  has  been 
utilized  on  this  configuration.  A review  of  various  cavitation  modeling  approaches  will 
then  be  presented.  Various  suggestions  and  modeling  ideas  will  then  be  presented  for 
analyzing  cavitation.  The  talk  concludes  with  a brief  summary  and  future  plans. 


14.6  Degree  Inducer-Impeller  with  Splitters 

This  view  graph  shows  the  computational  mesh  for  the  inducer-impeller  combination 
of  the  hub  surface.  The  inducer  and  impeller  rotate  at  the  same  RPM,  with  no  clocking 
in  between  them.  The  configuration  consists  of  a full  inducer  and  impeller  coupling  with 
a continuous  main  blade  that  proceeds  from  the  inducer  leading  edge  to  impeller  trailing 
edge.  A single  set  of  splitters  is  apparent  in  the  impeller.  The  figure  shows  the  blades 
and  splitters  projecting  out  of  the  hub  surface.  Although  there  are  four  total  main  blade 
passages,  only  one  was  considered  for  analysis  and  periodic  boundary  conditions  were 
applied  before  the  leading  edge  of  the  inducer  and  after  the  trailing  edge  of  the  impeller 
main  blade.  There  are  18  computational  nodes  in  between  the  main  blade  passages. 
There  are  about  110  points  along  the  flow  direction. 
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Summary  of  Pressure  Distribution  Over  the  Full  Configuration 

This  shows  the  pressure  distribution  on  the  full  configuration  (computational  results 
produced  with  a multiplicity  of  four)  with  splitters.  The  top  two  figures  show  the 
pressure  distribution  on  the  suction  and  pressure  sides  of  the  blades  and  splitters  and 
the  hub  surface  is  shown  shaded.  The  bottom  figures  shows  the  pressure  distributions 
on  the  hub,  centerline  and  shroud  surfaces  whereas  the  blades  and  splitters  shown  in 
the  grey  shade.  These  figures  indicate  the  existence  of  negative  pressure  near  tip  of  the 
suction  side  of  the  inducer  blade.  The  negative  pressures  are  due  to  the  fact  that  the 
CFD  model  assumes  a single-phase  flow  and  is  incapable  of  modeling  vapor  where  the 
local  static  pressure  becomes  lower  than  the  vapor  pressure. 


Current  Cavitation  Model  Applied  to  ATD  Configuration 

This  viewgraph  presents  an  approach  that  we  have  utilized  at  Pratt  Sz  Whitney  for 
getting  an  estimate  for  the  change  in  loading  attributed  to  cavitation.  The  first  step  is 
to  carry  out  a single  phase  3-D  CFD  solution  for  the  configuration  at  hand.  For  the  ATD 
configuration,  as  shown  in  the  previous  figure,  the  CFD  analysis  shows  that  the  static 
pressure  on  the  suction  side  near  the  tip  region  becomes  lower  than  local  vapor  pressure. 
Hence  the  effect  of  cavitation  must  only  manifest  itself  near  the  suction  side  tip  region 
and  gradually  diminish  as  one  move  towards  the  hub.  The  next  step  is  to  compute  a 
2-D  potential  flow  cavitation  sheet  model  corresponding  to  the  conditions  relevant  at 
the  blade  tip.  This  potential  flow  cavitation  sheet  model  has  been  developed  by  Penn 
State  University  for  P&W.  The  model  yields  cavitation  correction  factors,  Apsheet,  for 
pressure  loads  as  a function  of  tip  chord  length.  These  corrections  must  be  multiplied 
with  non-cavitated  blade  delta-pressures  to  obtain  cavitated  delta-pressures.  Note  that 
p is  used  here  for  pressure  corrections,  whereas  P is  used  for  static  pressures.  Next 
the  tip  sheet  cavitation  correction  factors  are  scaled  linearly  such  that  no  cavitation 
correction  is  needed  at  the  hub.  This  yields  a 2-D  distribution  of  the  correction  factors, 
ApCor,  on  the  blade  surface.  Thus  in 


A pcor  = K(s,r)Apsheet 

K(s,r)‘  represents  the  linear  scaling  transformation  as  a function  of  local  arc  length  s 
measured  from  the  inducer  LE  for  a streamline  at  a local  radius  r and  A pcor  represents 
the  multiplier  for  delta-pressure  at  all  blade  locations.  Once  the  functional  form  of 
the  delta-pressure  corrections  is  known,  it  is  multiplied  at  all  the  CFD  blade  locations 
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to  obtain  a new  distribution  A PCOT  for  blade  loads.  This  distribution  is  then  used  to 
compute  the  new  suction  side  pressure  corresponding  to  the  CFD  solution.  The  suction 
side  pressure  is  given  by: 

Ps  = PP  - A Pcor 

where  Pv  represents  the  pressure  on  the  pressure  side  of  the  blade.  As  a sanity  check,  one 
should  verify  that  the  corrected  pressure  on  the  suction  side  remains  above  the  vapor 
pressure.  If  this/ not  achieved,  then  the  2-D  sheet  model  should  be  repeated  for  other 
streamlines  and  the  cavitation  corrections  determined  by  a bilinear  interpolation. 


Distribution  of  Cavitation  Correction  Factors 

This  viewgraph  shows  the  functional  form  of  the  cavitation  correction  on  the  blade  tip 
as  a function  of  local  arc  length  along  the  blade.  The  effect  of  factors  greater  than  unity 
is  to  increase  the  pressure  load  downstream  of  the  bubble,  whereas  for  those  locations 
where  the  bubble  actually  exists,  the  effect  is  to  decrease  the  pressure  loads.  This  can  be 
thought  of  as  a blockage  effect  which  causes  the  variations  in  pressure.  The  figure  also 
shows  the  effect  of  linearly  scaling  the  correction  factors  and  its  diminishing  effect  as  one 
goes  from  tip  to  hub.  Note  that  this  procedure  should  be  regarded  as  rough  estimate 
of  the  changes  associated  with  cavitation.  For  this  reason  a more  elaborate  non-linear 
transfer  function  is  not  warranted.  Ideally  the  cavitation  correction  procedure  should  a 
part  of  the  CFD  numerics. 


Inducer  Blade  Delta  Pressures 

This  figure  shows  a projected  view  of  the  inducer  blade  in  the  x-r  plane.  The  top 
figure  shows  the  blade  pressure  loads  as  predicted  by  the  CFD  code  and  before  the 
cavitation  correction  is  applied.  The  cavitation  bubble  is  limited  to  the  red  region  near 
the  blade  tip.  The  lower  figure  shows  the  pressure  loads  after  the  cavitation  correction  is 
applied.  The  effect  of  this  correction  is  to  locally  decrease  the  load  due  to  the  blockage 
and  to  increase  the  loads  downstream  of  it.  The  saw-tooth  behavior  relates  to  the  axial 
coarseness  of  the  mesh  when  the  tip  correction  factors  are  linearly  transferred  to  other 
locations. 


Inducer  Blade  Suction-Side  Pressures  The  top  figure  shows  the  blade  suction 
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side  pressures  as  predicted  by  the  CFD  code  and  before  the  cavitation  correction  is 
applied  on  a projected  x-r  plane.  Note  that  the  pressures  are  substantially  negative 
in  the  cavity  region.  The  lower  figure  shows  the  suction  side  static  pressure  after  the 
cavitation  correction  is  applied.  In  this  case  all  the  pressures  have  become  positive. 
The  blockage  effect  is  felt  to  about  mid-span  location  at  the  tip.  The  blockage  effect 
diminishes  as  one  moves  from  tip  to  hub. 


Linear  Cavitation  Model 

I will  now  present  some  sheet  cavitation  models  that  one  can  consider.  Some  of  these 
have  been  considered  by  Penn  State.  Additional  ideas  are  presented  that  can  provide 
a more  accurate  and  simpler  treatment  of  the  cavitation  problem.  This  viewgraph 
describes  the  linear  cavitation  sheet  model.  In  such  a model  the  cavity  interface  is  treated 
as  a streamline  with  a constant  pressure  equal  to  the  vapor  pressure.  The  approach  is 
analogous  to  classical  linear  theory  used  in  the  potential  flow  codes.  Interface  conditions 
are  transferred  to  the  solid  body  and  zero  normal  velocity  condition  is  relaxed  at  the 
solid  boundary.  Thus  the  pressure  of  the  solid  boundary  adjacent  to  the  bubble  will  be 
equal  to  the  vapor  pressure.  The  relaxing  of  zero  normal  velocity  condition  at  the  solid 
boundary  simply  discards  all  flow  quantities  interior  to  the  cavity  and  correct  blockage 
effects  are  simulated  outside  of  the  cavity.  The  disadvantage  of  the  approach  is  that  it 
does  not  explicitly  treat  the  vapor  phase  and  that  the  approach  is  not  robust.  Since  the 
approach  handles  only  the  liquid  phase,  the  modeling  can  not  be  accurate  for  cryogenic 
fluids  operating  near  the  critical  point.  There  can  be  certain  cavity  termination  problems 
in  2-D  which  can  only  become  worse  in  three  spatial  dimensions. 


Non-Linear  Cavitation  Model 

In  the  non-linear  cavitation  model,  the  cavity  interface  is  again  regarded  as  a stream- 
line with  a specified  cavitation  pressure.  The  basic  difference  between  the  linear  and 
non-linear  models  is  that  the  computational  domain  evolves  with  the  solution  procedure 
as  a sequence  of  linear  solutions  followed  by  modifications  of  the  domain  to  accommo- 
date the  cavity  geometry.  Thus  the  computational  domain  is  regridded  after  locating  the 
positions  where  the  cross-over  below  the  vapor  pressure  occurs.  These  cross-over  points 
are  recomputed  subsequently  for  a better  definition  of  the  cavity.  One  can,  in  principle, 
iterate  to  a “correct”  cavity  geometry  via  these  multiple  passes.  The  disadvantage  of  the 


406 


approach  is  that  mesh  movement  may  introduce  significant  grid  distortions  and  the  grid 
lines  may  actually  cross  for  complicated  configurations.  Significant  man  power  efforts 
are  already  spend  in  gridding  complicated  configurations,  it  will  make  sense  to  retain 
these  meshes  for  a cavitating  flow.  Hence  a non-linear  cavitation  approach  that  holds 
the  base  grid  will  be  highly  desirable.  It  should  also  be  pointed  out  that  regridding  may 
only  be  suitable  for  box-like  domains  and  may  have  significant  problems  for  complex  3-D 
configurations.  For  example,  there  can  be  multiple  cavities  in  a domain,  or  a, cavity  may 
be  completely  embedded  in  the  domain  and  not  hooked  to  a physical  boundary.  There 
are  other  disadvantages  with  respect  to  the  physics  of  cavitation  sheet  models.  If  p > p„, 
then  the  fluid  remains  a liquid;  otherwise  the  fluid  is  in  the  vapor  phase.  No  dynamics  of 
the  vapor  phase  are  included  in  the  computations,  so  for  example,  the  change  in  volume 
in  flashing  to  a vapor,  or  the  collapse  in  volume  due  to  condensation  is  not  taken  into 
account.  These  effects  can  be  accounted  by  considering  a multi-phase  approach  with 
more  elaborate  description  of  physics.  However,  the  incorporation  of  two-phase  flows 
can  only  be  accommodated  in  a regridding  methodology  with  some  difficulty. 


Proposed  Cavitation  Model 

I will  now  present  a non-linear  two-phase  sheet  model  in  which  a structured  grid  does 
not  have  to  be  regenerated.  Firstly,  for  the  sake  of  all  subsequent  discussion,  consider 
the  static  pressure  to  be  referenced  from  vapor  pressure.  Thus  pressures  below  the  vapor 
pressure  are  regarded  as  negative  pressure  for  the  liquid.  First  carry  out  a single  phase 
(liquid)  CFD  solution  to  convergence.  Regard  cells  with  positive  pressure  on  all  nodes 
as  liquid  cells,  regard  cells  with  all  negative  node  pressures  as  vapor  cells  for  subsequent 
iterations,  and  regard  all  remaining  cells  as  interface  cells  where  both  liquid  and  vapor 
fluxes  will  be  carried  out  simultaneously.  The  third  step  is  to  determine  the  liquid  vapor 
interface  as  zero  crossing  for  pressure  on  the  edges  of  interface  cells.  Consider  a 2-D  cell 
acef  shown  in  the  viewgraph  for  reference.  Treat  the  liquid-vapor  interface  bd  locally 
as  a streamline.  This  means  that  one  should  discard  the  velocity  components  normal  to 
the  streamline.  Suppose  upb  is  the  predicted  value  of  velocity  at  the  zero  crossing  point 
b based  upon  linear  interpolation  of  the  edge  values  ac.  Let  ucb  be  the  corrected  value 
of  the  velocity  that  discards  the  velocity  component  normal  to  the  line  bd.  Then  redo 
the  flux  balance  for  the  liquid  part  of  the  cell  on  nodes  abdef.  Except  for  the  pressure 
terms  in  the  momentum  and  energy  equations,  the  interface  bd  yields  a zero  flux.  Thus 
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the  flux  balance  for  x-momentum  equation  on  the  interface  bd  is  simply: 

/ (pu2  + p)dy  - puvdx  = pvAyw. 

Jbd 

Based  upon  the  new  flux  balance,  one  can  compute  the  “distribution  formulas”  for 
nodes  aef.  These  distribution  relations  yield  contribution  of  the  cell  flux  balance  to 
the  corresponding  nodes  (See  Ni’s  classical  AIAA  Journal  Paper,  1985).  A similar  flux 
balance  can  be  carried  out  for  the  vapor  part  of  the  cell.  One  can,  in  fact,  regard  the 
interface  cell  nodes  to  be  multiply  defined  for  vapor  and  liquid  parts.  Thus  the  liquid 
part  will  have  real  fluid  values  at  nodes  ae f and  fictitious  values  at  node  c.  The  flux 
balance  can  be  carried  out  independently  based  upon  the  real  and  fictitious  values  and 
the  nodes  updated  separately.  This  approach  constitutes  a simple  two-phase  approach 
to  cavitation  without  the  mesh  regridding  option.  Note  that  I have  not  specifically 
stated  anything  about  the  physics  of  two-phase  flows  yet.  A more  elaborate  two-phase 
approach  without  appropriate  physics  may  only  yield  qualitatively  correct  answers.  This 
follows  on  the  next  viewgraph. 


Additional  Modeling  Enhancements 

A more  elaborate  equation  of  state  which  is  approximately  valid  in  the  gas/liquid 
two-phase  region,  as  well  as  at  moderately  high  temperatures  is  needed.  A suitable 
model  that  describes  the  law  of  corresponding  states  can  be  easily  implemented  for  the 
thermal  equation  of  state  (that  relates  pressure,  density  and  temperature).  This  has 
to  be  coupled  with  a more  representative  caloric  equation  of  state  (that  relates  specific 
heat  to  temperature).  One  has  to  also  employ  a more  accurate  energy  equation  to 
describe  two-phase  flow  physics.  This  basically  includes  latent  heat  of  vaporization  for 
the  vapor  phase.  Certain  cryogenic  fluids,  with  critical  temperatures  below  50  K,  exhibit 
quantum  effects,  and  do  not  necessarily  satisfy  the  principle  of  corresponding  states.  A 
more  elaborate  equation  of  state,  possibly  in  the  form  of  look  up  tables,  may  be  needed 
for  these  fluids.  Note  that  both  liquid  hydrogen  and  helium  exhibit  these  quantum 
effects.  Liquid  hydrogen  is  a common  fluid  for  rocket  engine  fuel  pumps.  Lastly,  one 
can  utilize  mesh  adaptation  to  locally  subdivide  the  cells  in  the  vicinity  of  liquid- vapor 
interface  to  impart  additional  spatial  resolution.  This  will  eliminate  the  need  to  carry 
out  multi-valued  flux  balance,  providing  that  the  embedded  cells  are  reasonably  resolved. 
Note  that  embedded  mesh  adaptation  does  not  regrid  the  “base”  mesh,  the  skewness  of 
the  embedded  meshes  will  only  be  limited  by  the  skewness  of  the  coarse  mesh,  and  that 
cells  can  never  cross  if  the  base  mesh  is  not  crossed  over.  Pratt  &:  Whitney  already  has 
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suitable  data  base  structure  for  handling  this  approach. 


Closure 

A CFD  solution  for  ATD  inducer /impeller  pump  configuration,  with  splitters,  has 
been  presented  at  the  design  point.  This  solution  points  to  the  existence  of  a cavitation 
bubble  in  the  domain.  Current  cavitation  approach  employs  a separate  potential  flow 
cavitation  sheet  model  that  yields  cavitation  correction  factors  for  the  CFD  predicted 
pressure  loads.  Penn  State  University  is  currently  studying  linear  and  non-linear  ap- 
proaches for  handling  3-D  configurations.  If  a suitable  linear  cavitation  approach  can 
constitute  a simple  and  effective  model,  we  will  implement  such  a model  in  the  Pratt  h 
Whitney  design  code.  Additional  enhancements  for  two-phase  flows  and  more  accurate 
physics  will  be  considered  as  part  of  a long-term  effort. 
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AN  INDUCER  CFD  SOLUTION  AND  EFFECTS 
ASSOCIATED  WITH  CAVITATION 
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West  Palm  Beach,  Florida 
April  20, 1993 

Presented  At  Workshop  for  CFD 
Applications  in  Rocket  Propulsion,  NASA— MSFC 


□ Atd  inducer/impeller  CFD  Analysis 
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ORIGINAL  PAGE  IS 
OF  POOR  QUALITY 


CURRENT  CAVITATION  PROCEDURE  APPLIED 
TO  ATD  CONFIGURATION 
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Q As  a sanity  check  verify  that  the  corrected  suction  pressure 
remains  above  the  vapor  pressure. 


ATD  14.6  I n d u c e r- I m p e 1 1 e r CFD  Predictions 
Spanwise  Variation  of  Cavitation  Correction  Factor 
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ATD  14.6  Inducer-Impeller  CFD  Predictions 
Inducer  Blade  Delta  Pressures 
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Corrected  for  Cavitation 


ATD  14.6  Inducer-Impeller  CFD  Predictions 
Inducer  Blade  Suction-Side  Pressure 
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Corrected  for  Cavitation 
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-LINEAR  CAVITATION  MODEL 
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PROPOSED  CAVITATION  MODEL 
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□ Liquid  and  Vapor  Nodes  of  Interface  Cells  Will  Vary  Independently. 

□ This  Constitutes  a Non-linear  Multi -phase  approximation  to 
Cavitation  Without  the  Mesh  Regridding  Option. 


ADDITIONAL  MODELING  ENHANCEMENTS 
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out  multi-valued  flux  balance. 


CLOSURE 
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□ Additional  enhancements  will  be  considered  in  the  future. 


N95- 23632 


CURRENT  STATUS  IN  CAVITATION  MODELING 


Ashok  K.  Singhal  and  Ram  K.  Awa 
CFD  Research  Corporation,  Huntsville,  AL 


for  presentation  at 

Open  Forum  on  "Turbomachinery-Pumps-Cavitation" 


11th  Workshop  for  CFD  Applications  in  Rocket  Propulsion 
NASA  Marshall  Space  Flight  Center,  AL 


April  20-22, 1993 

Cavitation  is  a common  problem  for  many  engineering  devices  in  which  the  main 
working  fluid  is  in  liquid  state.  In  turbomachinery  applications,  cavitation  generally 
occurs  on  the  inlet  side  of  pumps.  The  deleterious  effects  of  cavitation  include: 
lowered  performance,  load  asymmetry,  erosion  and  pitting  of  blade  surfaces, 
vibration  and  noise,  and  reduction  of  the  overall  machine  life. 


Cavitation  models  in  use  today  range  from  rather  crude  approximations  to 
sophisticated  bubble  dynamics  models.  Details  about  bubble  inception,  growth  and 
collapse  are  relevant  to  the  prediction  of  blade  erosion,  but  are  not  necessary  to 
predict  the  performance  of  pumps.  An  engineering  model  of  cavitation  is  proposed 
to  predict  the  extent  of  cavitation  and  performance.  The  vapor  volume  fraction  is 
used  as  an  indicator  variable  to  quantify  cavitation.  A two-phase  flow  approach  is 
employed  with  the  assumption  of  the  thermal  equilibrium  between  liquid  and 
vapor.  At  present  velocity  slip  between  the  two  phases  is  selected.  Preliminary 
analyses  of  2D  flows  shows  qualitatively  correct  results. 
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PA-93-8PR/01 


BOILING  VS  CAVITATION 
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due  to  lowering  of 
hydrodynamic  pressure 
under  adiabatic  conditions. 


HIERARCHY  OF  CAVITATION  MODELS 
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Selected  Approach:  Level  3 


DENSITY  CHANGES  IN  ENGINEERING  FLOWS 


Smaller  sound  speed  slows  down  information  propagation 
Cavitating  flows  offer  a very  stringent  challenge  to  numerical  schemes 


TEST  PROBLEM 
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Study  2:  Effect  of  geometry  change  (rounded 

entrance;  ylR  = 1 5%) 


CATERPILLAR  INJECTOR  ORIFICE 
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c.  VAPOR  FRACTION 


(DHOOOOOOOOH 
H rl  rl  I I I I I I I I 

OOO®®®^®®®41 

OOOrHC©IXl<Ni5>^<oy5 

OOOOr^OLfSfOiNHfD 


c^vofo  r-.  rp  (\j  o yD  (N 

rsjfO'TLDLnor-r-coO 

C^tOr-^Lfl'COfMHO 



©ooooooooo 


Hi  S 

o m 

— II 

U. 

E ® 

O .E 

CO  0- 


QC  < I 
< 0)  ■ 

j 2 I 

E w I 

oc  > If 

Li!  U 

h-  *0  I 

< S 1 

O 00  ft 


oooooooooo 

+ + + + + + + + + + 

^5TrOfOC©CNr^r4'J,0>lCO 


fO  (N  <N  (N 


H H H CTi  O 


c.  VAPOR  FRACTION  ■ l . 4 9 e - o 8 


CATERPILLAR  INJECTOR  ORIFICE 
Steady  State  Analysis  Pinlet  = 10.0  bar 


c.  VAPOR  FRACTION 


CATERPILLAR  INJECTOR  ORIFICE 

Effect  of  Edge  Shape  on  Void  Fraction  Pinj=4  bar 
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Rounded  Edge 
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3)  Validation  under  a variety  of  flow  conditions. 


N95- 23633 


A GENERALIZED  EULERIAN-LAGRANGIAN  ANALYSIS,  WITH  APPLICATION  TO 
LIQUID  FLOWS  WITH  VAPOR  BUBBLES^ 


Frederik  J.  de  Jong  and  Meyya  Meyyappan 

Scientific  Research  Associates,  Inc:  j / 

Glastonbury,  CT 


Presented  at  the 

Workshop  for  Computational  Fluid  Dynamics  Applications 
in  Rocket  Propulsion 

April  20-22, 1993 


ABSTRACT 


Under  a NASA  MSFC  SBIR  Phase  II  effort  an  analysis  has  been  developed  for  liquid  flows  with 
vapor  bubbles  such  as  those  in  liquid  rocket  engine  components.  The  analysis  is  based  on  a 
combined  Eulerian-Lagrangian  technique,  in  which  Eulerian  conservation  equations  are  solved  for  the 
liquid  phase,  while  Lagrangian  equations  of  motion  are  integrated  in  computational  coordinates  for 
the  vapor  phase.  The  novel  aspect  of  the  Lagrangian  analysis  developed  under  this  effort  is  that  it 
combines  features  of  the  so-called  particle  distribution  approach  with  those  of  the  so-called  particle 
trajectory  approach  and  can,  in  fact,  be  considered  as  a generalization  of  both  of  those  traditional 
methods.  The  result  of  this  generalization  is  a reduction  in  CPU  time  and  memory  requirements. 
Particle  time  step  (stability)  limitations  have  been  eliminated  by  semi-implicit  integration  of  the  particle 
equations  of  motion  (and,  for  certain  applications,  the  particle  temperature  equation),  although 
practical  limitations  remain  in  effect  for  reasons  of  accuracy.  The  analysis  has  been  applied  to  the 
simulation  of  cavitating  flow  through  a single-bladed  section  of  a labyrinth  seal.  Models  for  the 
simulation  of  bubble  formation  and  growth  have  been  included,  as  well  as  models  for  bubble  drag 
and  heat  transfer.  The  results  indicate  that  bubble  formation  is  more  or  less  "explosive*:  for  a given 
flow  field,  the  number  density  of  bubble  nucleation  sites  is  very  sensitive  to  the  vapor  properties  and 
the  surface  tension.  The  bubble  motion,  on  the  other  hand,  is  much  less  sensitive  to  these 
properties,  but  is  affected  strongly  by  the  local  pressure  gradients  in  the  flow  field.  In  situations  where 
either  the  material  properties  or  the  flow  field  are  not  known  with  sufficient  accuracy,  parametric 
studies  can  be  carried  out  rapidly  to  assess  the  effect  of  the  important  variables.  Future  work  will 
include  application  of  the  analysis  to  cavitation  in  inducer  flow  fields. 


This  work  was  supported  by  NASA  Marshail  Space  Flight  Center  under  Contract  NAS8-38959. 
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PRECEDING, 

PAGE.  7-^ 


. 6LA.NK  NOT  FILMED 
INTENTIONALLY  BLANK 


OBJECTIVE: 

DEVELOPMENT  OF  AN  ANALYSIS  FOR  LIQUID  FLOWS  WITH  VAPOR 
BUBBLES  (SUCH  AS  THOSE  IN  BEARINGS,  SEALS,  AND  PUMPS)  FOR 
USE  IN  DESIGN  OF  LIQUID  ROCKET  ENGINE  COMPONENTS 


APPROACH: 

USE  A COMBINED  EULER1AN-LAGRANGIAN  ANALYSIS 

- CONTINUOUS  (LIQUID)  PHASE  TREATED  BY  SOLVING  EULERIAN 
CONSERVATION  EQS.  (N-S  EQS.)  WITH  DISCRETE  PHASE  SOURCE  TERMS 

- DISCRETE  (VAPOR  BUBBLE)  PHASE  TREATED  BY  INTEGRATING 
LAGRANGIAN  EQUATIONS  OF  MOTION  IN  COMPUTATIONAL  COORDINATES 


LAGRANGIAN  ANALYSIS  - METHODOLOGY 


GENERALIZATION  OF  PARTICLE  TRAJECTORY  AND  PARTICLE 
DISTRIBUTION  MODELS 

=*  REDUCED  CPU  TIME  AND  STORAGE  REQUIREMENTS 

STABLE  INTEGRATION  OF  LAGRANGIAN  EQUATIONS  OF 
MOTION  IN  COMPUTATIONAL  COORDINATES 


Scientific 

Research 

Associates 


NO  PARTICLE  TIME  STEP  LIMITATION  OTHER 
THAN  FOR  ACCURACY 
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LAGRANGIAN  ANALYSIS  - FEATURES 


LAGRANGIAN  DISCRETE  PHASE  ANALYSIS  INTERFACES 
WITH  EULERIAN  CONTINUOUS  PHASE  ANALYSIS 
VIA  ONE  SINGLE  SUBROUTINE 

GENERALIZED  ARRAY  ADDRESSING  AND  INDEX  / VARIABLE 
NUMBERING 

=>  LAGRANGIAN  ANALYSIS  CAN  BE  HOOKED  UP  TO 
A VARIETY  OF  EULERIAN  CODES  ("FLOW  SOLVERS”) 

CODE  IS  MODULAR  W.R.T.  PHYSICAL  MODELS 


APPLICATION  TO  VAPOR  BUBBLES 


MODELS  NEEDED  FOR  BUBBLE  FORMATION, 
GROWTH  / COLLAPSE  AND  MOTION 


CRITICAL  PHYSICAL  PARAMETERS: 

- SURFACE  TENSION 

- VAPOR  PRESSURE  - TEMPERATURE  RELATION 
(CLAUSIUS-CLAPEYRON  EQUATION) 


Scientific 

Research 

Associates 
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BUBBLE  FORMATION 

• NUCLEATION  DUE  TO  LOCAL  P DROP; 

USE  HOMOGENEOUS  NUCLEATION  MODEL 
(KATZ  AND  BLANDER) 


z NUCLEI  x _ / 2a  x 1/2  r 
( IX  3q  )= n (tit)  “P[- 


•16Ka3 

3kT(P-  Pl)' 


NUMBER  OF  NUCLEATION  SITES  PER  UNIT  VOLUME 
AND  TIME 

EQUILIBRIUM  VAPOR  PRESSURE 

PRESSURE  IN  LIQUID 

TEMPERATURE  OF  THE  LIQUID 

SURFACE  TENSION 

BOLTZMANN  CONSTANT 

MOLECULAR  WEIGHT 

MOLECULAR  NUMBER  DENSITY  IN  LIQUID 


BUBBLE  GROWTH 


AFTER  THE  INITIAL  GROWTH  PERIOD  (WHICH  IS  INERTIA 
CONTROLLED)  THE  PRESSURE  INSIDE  THE  BUBBLE  IS 
DETERMINED  BY  THE  YOUNG-LAPLACE  EQUATION 


PB  - P L = 2a  / R 


THE  BUBBLE  TEMPERATURE  IS  ASSUMED  TO  BE 

EQUAL  TO  THE  SATURATION  TEMPERATURE  AT  THE  PRESSURE  £ 


BUBBLE  GROWTH  IS  DETERMINED  BY  THE  ENERGY  EQUATION 
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ENERGY  EQUATION  FOR  BUBBLE 


q = m cpTp  + m h,g 

cp  : VAPOR-PHASE  SPECIFIC  HEAT 
hjg  : LATENT  HEAT  OF  EVAPORATION 
m : BUBBLE  MASS  GROWTH  RATE  ("EVAPORATION  RATE”) 
T*  : RATE  OF  CHANGE  OF  BUBBLE  TEMPERATURE 


* / \ -* 

tp=(^Lu*vp 


HEAT  TRANSFER  RATE  TO  BUBBLE 

q s Nu  D k (T  - Tp) 

Nu 

: NUSSELT  NUMBER 

D 

: BUBBLE  DIAMETER 

K 

: THERMAL  CONDUCTIVITY  OF  THE  LIQUID 

T 

: LIQUID  TEMPERATURE 

TP 

: BUBBLE  TEMPERATURE 

Scientific 

Research 

vrr 
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BUBBLE  EQUATION  OF  MOTION 


p : LIQUID  DENSITY 
p : VAPOR  DENSITY 

p 

Aa  : ADDED  MASS  COEFFICIENT  (=1  FOR  STOKES  FLOW) 
m : BUBBLE  MASS 
2 : BUBBLE  POSITION  VECTOR 

F : FORCE  ON  BUBBLE 

- DRAG  FORCE 

- BUOYANCY 

- PRESSURE  GRADIENT 


DRAG  COEFFICIENT 


FOR  BUBBLES 

co(=_Tt)  =14.9 /Re 0,78  Re >2 
npazuz 


= 16  /Re 


Re  <2 


GOOD  FIT  TO  NUMERICAL  SOLUTION 

COMPARES  WELL  WITH  EXPERIMENTS,  WIDELY  USED 

FOR  SOLID  SPHERES 

C„=  ■—  [1  +0.15  Re0'687]  Re  <1000 
Re 

= 0.438,  Re  >1000 

NOTE:  FOR  SAME  a AND  u,  DRAG  ON  A BUBBLE 
IS  SMALLER  (FACTOR  6 FOR  Re  = 900) 


Scientific 

Research 
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COMPUTATIONAL  TEST  CASE 


SINGLE  BLADE  LABYRINTH  SEAL 

APPROXIMATES  ONE  BLADE  OF  THE  LOX  PREBURNER 
BOOST  PUMP  IMPELLER  SEAL  IN  THE  SSME 


FLOW  CONDITIONS: 

INLET  PRESSURE  : 3.67  MPa  (532  psi) 
INLET  TEMPERATURE  : 130  K (234  °R) 


DENSITY 

PRESSURE  DROP 


1088  kg/m3  (1 67.915/1^ ) 
1.47  MPa  (213  psi) 
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GEOMETRY 


J* V 70.0 
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SEAL  GAP:  1.524  x 10  m (0.006  Inch) 
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COMPUTATIONAL  GRID 


AXIAL  VELOCITY 


BUBBLE  TRAJECTORIES 


BUBBLE  SIZE 
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N95-  23634 


On  the  Accuracy  of  CFD-Based  Pressure  Drop  r[„  f $ 
Predictions  for  Right- Angle  Ducts 


A.  Brankovic 
Pratt  & Whitney,  Florida 


The  predictive  capability  of  computational  fluid  dynamics  (CFD)  codes  for  turbu- 
lent flow  through  curved  ducts  is  of  significant  importance  to  the  design  and  perfor- 
mance analysis  of  modern  rocket  engine  flowpaths.  Code  calibration  and  validation 
studies  for  this  class  of  flow  are  desireable  to  estimate  the  performance  margin  and 
operating  range  of  components  designed  using  Navier-Stokes  methods.  Parametric 
experimental  studies  such  as  that  of  Weske  (NACA  ARR  W-39)  provided  a wealth  of 
performance  data  for  the  design  of  single-  and  compound  elbow  configurations  with 
various  cross-sections,  curvature  and  aspect  ratios  at  varying  Reynolds  numbers.  In 
that  work,  the  majority  of  data  is  presented  in  the  form  of  loss  coefficients,  character- 
izing pressure  losses  due  to  duct  curvature,  and  including  losses  due  to  wall  friction. 
Using  measured  friction  coefficients,  losses  of  equivalent  straight  lengths  of  duct  are 
subtracted,  resulting  in  performance  curves  useful  for  design  computations.  These 
data  are  currently  used  in  a CFD-based  parametric  study  covering  a broad  range 
of  operating  conditions.  Of  particular  interest  for  the  accuracy  of  CFD  predictions 
are  the  effects  on  pressure  loss  due  to  inlet  boundary  layer  thickness  (dependent  on 
upstream  development  length),  and  the  wall  treatment  for  the  turbulence  equations 
(conventional  wall  functions  vs.  wall  integration  using  a two- layer  model).  The  ex- 
perimental data  are  reassessed  in  the  form  of  an  error  analysis,  and  are  compared 
with  CFD  predictions  for  18  computational  cases.  Grid-independence,  grid  spacing, 
and  convergence  requirements  of  the  cases  are  discussed.  Conclusions  regarding  the 
relative  importance  of  the  parametric  variables  will  be  presented. 
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ON  THE  ACCURACY  OF  CFD -BASED 
PRESSURE  DROP  PREDICTIONS  FOR 
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EXPERIMENTAL  DATASETS 


Reference:  Weske,  J.R.,  ’’Pressure  Loss  in  Ducts  with 

Compound  Elbows”,  NACA  Wartime  Report  W— 39 
February  1943. 


EXPERIMENTAL  GEOMETRY 
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CFD  NUMERICS 

Grid  Independence 
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Similar  results  obtained  for  wall  function  and  wall  integration  models 


PREDICTED  vs.  MEASURED  RESULTS 

Influence  of  Duct  Aspect  Ratio  at  If  = 10 5 
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Data  shows  distinct  minimum  when  plotted  as 


COMPUTED  FLOW  STRUCTURE 

Influence  of  Duct  Aspect  Ratio  at  Re  = 105 
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Shown  at  48  1 Hownshpnm  of  plhow  exit  plane 


PREDICTED  vs.  MEASURED  RESULTS 

Influence  of  Duct  Aspect  Ratio  at  Re  = 6 xlO5 
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Minimum  in  data  not  as  pronounced  when  plotted  as 


PREDICTED  vs.  MEASURED  RESULTS 

Influence  of  Duct  Curvature  Ratio  at  Re  = 6 xlO5 
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Dnct  Cinvature  Ratio 


FLOW  SEPARATION  DETAILS 

Turbulence  Model  Influences  Shape,  Location  Of  Separation  Bubble 
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N95- 23635 


CFD  MODELING  OF  TURBULENT  DUCT  FLOWS  FOR 
COOLANT  CHANNEL  ANALYSIS 

Ronald  J.  Ungewitter  and  Daniel  C.  Chan 
Rocketdyne  Div.  /Rockwell  International 
Canoga  Park,  Ca  91304 

ABSTRACT  j _ (J 

The  design  of  modern  liquid  rocket  engines  requires  the  analysis  of  chamber 
coolant  channels  to  maximize  the  heat  transfer  while  minimizing  the  coolant  flow.  Coolant 
channels  often  do  not  remain  at  a constant  cross  section  or  at  uniform  curvature.  New 
designs  require  higher  aspect  ratio  coolant  channels  than  previously  used.  To  broaden  the 
analysis  capability  and  to  complement  standard  analysis  tools  an  investigation  on  the 
accuracy  of  CFD  predictions  for  coolant  channel  flow  has  been  initiated.  Validation  of 
CFD  capabilities  for  coolant  channel  analysis  will  enhance  the  capabilities  for  optimizing 
design  parameters  without  resorting  to  extensive  experimental  testing.  The  eventual  goal  is 
to  use  CFD  to  determine  the  flow  fields  of  unique  coolant  channel  designs  and  therefore 
determine  critical  heat  transfer  coefficients. 

In  this  presentation  the  accuracy  of  a particular  CFD  code  is  evaluated  for  turbulent 
flows.  The  first  part  of  the  presentation  is  a comparison  of  numerical  results  to  existing 
cold  flow  data  for  square  curved  ducts  (NASA  CR-3367,  "Measurements  of  Laminar  and 
Turbulent  Flow  in  a Curved  Duct  with  Thin  Inlet  Boundary  Layers").  The  results  of  this 
comparison  show  good  agreement  with  the  relatively  coarse  experimental  data.  The  second 
part  of  the  presentation  compares  two  cases  of  higher  aspect  ratio  channels  (AR=2.5,10)  to 
show  changes  in  axial  and  secondary  flow  strength.  These  cases  match  experimental  work 
presently  in  progress  and  will  be  used  for  future  validation.  The  comparison  shows 
increased  secondary  flow  strength  of  the  higher  aspect  ratio  case  due  to  the  change  in  radius 
of  curvature.  The  presentation  includes  a test  case  with  a heated  wall  to  demonstrate  the 
program's  capability.  The  presentation  concludes  with  an  outline  of  the  procedure  used  to 
validate  the  CFD  code  for  future  design  analysis. 
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FLOW  THROUGH  A 90  SQUARE  BENT  DUCT:  Re-40,000 
AXIAL  VELOCITY  MAGNITUDE:  Theta-60  deg 


REACT3D  CFD  RESULTS 
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FLOW  THROUGH  A 90  SQUARE  BENT  DUCT:  Re  = 40,000 
AXIAL  VELOCITY  MAGNITUDE:  Y=0.25 


REACT3D  CFD  RESULTS, 
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PRESSURE  VALUES  ARE  ALONG  THE  OUTSIDE  WALL,  Z=0.0,  AND 
SIDE  WALLS,  R*=0.1,  0.3, .05,  0.7,  0.9 
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REACT3D  SOLUTION  OF  60  DEG.  RECT.  DUCT 

ASPECT  RATIO=  10,  RE=  127,000,  THETA=0 
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REACT3D  SOLUTION  OF  60  DEG.  RECT.  DUCT 

ASPECT  RATIO  = 2.5,  RE=  100,000,  THETA  = 0 
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REACT3D  SOLUTION  OF  60  DEG.  RECT.  DUCT 

ASPECT  RATI0  = 2.5,  RE=  100,000,  THETA=60 
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REACT3D  SOLUTION  OF  60  DEG.  RECT.  DUCT 

ASPECT  RATIO  = 10,  EE=  127,000,  THETA=60 
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REACT3D  SOLUTION  OF  60  DEG.  RECT.  DUCT 

ASPECT  RATIO  = 2.5,  RE=  100,000,  THETA=60 
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REACT3D  SOLUTION  OF  60  DEG.  RECT.  DUCT 

ASPECT  RATIO=  10,  RE=  127,000,THETA=60 
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FLOW  AND  HEAT  TRANSFER  IN  180-DEGREE  TURN  SQUARE  DUCTS 
- EFFECTS  OF  TURNING  CONFIGURATION  AND  SYSTEM  ROTATION 

Ten-See  Wang 

Computational  Fluid  Dynamics  Branch 
NASA  - Marshall  Space  Flight  Center 
Marshall  Space  Flight  Center,  AL  35812 

Ming-King  Chyu 

Department  of  Mechanical  Engineering 

Carnegie  Mellon  University  Q /# 

Pittsburgh,  PA  15213  > 


ABSTRACT 

Forced  flow  through  channels  connected  by  sharp  bends  is 
frequently  encountered  in  various  rocket  and  gas  turbine  engines. 
For  example,  the  transfer  ducts,  the  coolant  channels  surround  the 
combustion  chamber,  the  internal  cooling  passage  in  a blade  or 
vane,  the  flow  path  in  the  fuel  element  of  a nuclear  rocket  engine, 
the  flow  around  a pressure  relieve  valve  piston,  and  the 
recirculated  base  flow  of  multiple  engine  clustered  nozzles. 
Transport  phenomena  involved  in  such  a flow  passage  are  complex  and 
considered  to  be  very  different  from  those  of  conventional  turning 
flow  with  relatively  mild  radii  of  curvature.  While  previous 
research  pertaining  to  this  subject  has  been  focused  primarily  on 
the  experimental  heat  transfer,  very  little  analytical  work  is 
directed  to  understanding  the  flowfield  and  energy  transport  in  the 
passage.  Therefore,  the  primary  goal  of  this  paper  is  to  benchmark 
the  predicted  wall  heat  fluxes  using  a state-of-the-art 
computational  fluid  dynamics  (CFD)  formulation  against  those  of 
measurement  for  a rectangular  turn  duct.  Other  secondary  goals 
include  studying  the  effects  of  turning  configurations,  e.g.,  the 
semi-circular  turn,  and  the  rounded-corner  turn,  and  the  effect  of 
system  rotation.  The  computed  heat  fluxes  for  the  rectangular  turn 
duct  compared  favorably  with  those  of  the  experimental  data.  The 
results  show  that  the  flow  pattern,  pressure  drop,  and  heat 
transfer  characteristics  are  different  among  the  three  turning 
configurations,  and  are  substantially  different  with  system 
rotation.  Also  demonstrated  in  this  work  is  that  the  present 
computational  approach  is  quite  effective  and  efficient  and  will 
be  suitable  for  flow  and  thermal  modeling  in  rocket  and  turbine 
engine  applications. 
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Flow  and  Heat  Transfer  in  180-Degree  Turn  Square  Ducts 
- Effects  of  Turning  Configuration  and  System  Rotation 
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■ To  examine  the  influence  of  turn  configuration  and  System 
rotation  on  the  overall  transport  phenomena  in  the  passage 
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The  evolution  of  the  secondary  flow  pattern  for  a 180-degree  turn  straight-corner  square  duct 
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The  evolution  of  the  secondary  flow  pattern  for  a 180-degree  turn  round-corner  square  duct 
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The  evolution  of  the  secondary  flow  pattern  for  a 180-degree  turn  circular-corner  square  duct 
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The  evolution  of  the  secondary  flow  pattern  for  a 180-degree  turn  straight-corner  square  duct 
with  rotation 
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The  averaged  total  Nusselt  numbers 
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Conclusions 
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• Circular-corner  turn  has  highest  post-turn  heat 
transfer 

• System  rotation  enhances  heat  transfer 


N95- 23637 


Methodology  for  CFD  Design  Analysis 
of  National  Launch  System 
Nozzle  Manifold 

/ } 

' ■ o / 


Scot  L.  Haire 

Pratt  & Whitney,  West  Palm  Beach,  FL 


ABSTRACT 


The  current  design  environment  dictates  that  high  technology  CFD 
(Computational  Fluid  Dynamics)  analysis  produce  quality  results  in  a 
timely  manner  if  it  is  to  be  integrated  into  the  design  process.  The 
design  methodology  outlined  describes  the  CFD  analysis  of  an  NLS  (Na- 
tional Launch  System)  nozzle  film  cooling  manifold.  The  objective  of 
the  analysis  was  to  obtain  a qualitative  estimate  for  the  flow  distribution 
within  the  manifold.  A complex,  3D,  multiple  zone,  structured  grid  was 
generated  from  a 3D  CAD  file  of  the  geometry.  An  Euler  solution  was 
computed  with  a fully  implicit  compressible  flow  solver.  Post  processing 
consisted  of  full  3D  color  graphics  and  mass  averaged  performance.  The 
result  was  a qualitative  CFD  solution  that  provided  the  design  team  with 
relevant  information  concerning  the  flow  distribution  in  and  performance 
characteristics  of  the  film  cooling  manifold  within  an  effective  time  frame. 
Also,  this  design  methodology  was  the  foundation  for  a quick  turnaround 
CFD  analysis  of  the  next  iteration  in  the  manifold  design. 
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West  Palm  Beach,  Florida 
April  20, 1993 

Presented  At  Workshop  for  CFD 
Applications  in  Rocket  Propulsion,  NASA— MSFC 


NLS  MANIFOLD  ANALYSIS 
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ANALYSIS  OBJECTIVE 
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STME  FILM/CONVECTIVE  COOLED  NOZZLE 
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ADVANCED  MULTI-PHASE  FLOW  CFD  MODEL  DEVELOPMENT  FOR 
SOLID  ROCKET  MOTOR  FLOWFIELD  ANALYSIS* 


Paul  Liaw,  Y.S.  Chen,  and  H.  M.  Shang 
Engineering  Sciences,  Inc. 
Huntsville,  AL 


D.  Doran 

ED  32,  NASA  Marshall  Space  Flight  Center 


Xb 


ABSTRACT 


It  is  known  that  the  simulations  of  solid  rocket  motor  internal  flow  field  with  AL-based  propellants 
require  complex  multi-phase  turbulent  flow  model.  The  objective  of  this  study  is  to  develop  an  advanced 
particulate  multi-phase  flow  model  which  includes  the  effects  of  particle  dynamics,  chemical  reaction  and 
hot  gas  flow  turbulence.  The  inclusion  of  particle  agglomeration,  particle/gas  reaction  and  mass  Iransfer, 
particle  collision,  coalescence  and  breakup  mechanisms  in  modeling  the  particle  dynamics  will  allow  the 
proposed  model  to  realistically  simulate  the  flowfield  inside  a solid  rocket  motor. 

The  Finite  Difference  Navier-Stokes  numerical  code  FDNS  is  used  to  simulate  the  steady-state 
multi-phase  particulate  flow  field  for  a 3-zone  2-D  axisymmetric  ASRM  model  and  a 6-zone  3-D  ASRM 
model  at  launch  conditions.  The  2-D  model  includes  aft-end  cavity  and  submerged  nozzle.  The  3-D  model 
represents  the  whole  ASRM  geometry,  including  additional  grain  port  area  in  the  gas  cavity  and  two 
inhibitors. 


FDNS  is  a pressure  based  finite  difference  Navier-Stokes  flow  solver  with  time-accurate  adaptive 
second-order  upwind  schemes,  standard  and  extended  k-e  models  with  compressibility  corrections,  multi- 
zone body-fitted  formulations,  and  turbulence/particle  interaction  model.  Eulerian/Lagrangian  multi-phase 
solution  method  is  applied  for  multi-zone  mesh.  To  simulate  the  chemical  reaction,  penalty  function 
corrected  efficient  finite-rate  chemistry  integration  method  is  used  in  FDNS.  For  the  AL  particle 
combustion  rate,  the  Hermsen  correlation  is  employed.  To  simulate  the  turbulent  dispersion  of  particles, 
the  Gaussian  probability  distribution  with  standard  deviation  equal  to  (2k/3)ic  is  used  for  the  random 
turbulent  velocity  components. 

The  flow  field  in  the  aft-end  cavity  of  the  ASRM  is  analyzed  to  investigate  its  significant  impact  on 
the  operation  of  the  motor  as  well  as  its  performance.  It  is  known  that  heat  flux  and  the  pressure 
distributions  in  this  region  will  cause  recirculation  and  influence  the  design  requirements.  Chemical 
reaction  of  gas  flow  is  a factor  affecting  the  performance  of  the  ASRM.  An  accurate  analysis  for 
chemically  reacting  flow  is  therefore  important  in  the  design  of  the  ASRM.  Twelve  gas  elements  (H20, 
02,  H2,  O,  H,  OH,  CO,  C02,  CL,  CL2,  HCL,  and  N2)  were  considered  for  the  chemical  reaction  in  present 
study.  For  multi-phase  calculations,  the  particulate  phase  was  injected  at  the  propellant  grain  surface. 
The  particulate  phase  was  assumed  to  be  aluminum  oxide  (A1203)  only.  The  mass  fraction  of  the 
particulate  phase  was  assumed  to  be  53%  of  the  mixture. 

The  computational  results  reveal  that  the  flow  field  near  the  juncture  of  aft-end  cavity  and  the 
submerged  nozzle  is  very  complex.  The  effects  of  the  turbulent  particles  affect  the  flow  field  significantly 
and  provide  better  prediction  of  the  ASRM  performance.  The  multi-phase  flow  analysis  using  the  FDNS 
code  in  the  present  research  can  be  used  as  a design  tool  for  solid  rocket  motor  applications. 

* This  work  is  supported  by  NASA  Marshall  Space  Flight  Center  under  Contract  NAS8-39398. 
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Abstract 


The  Advanced  Solid  Rocket  Motor  (ASRM)  program  was  initiated  by  NASA  in 
response  to  the  need  for  a new  generation  rocket  motor  capable  of  providing 
increased  thrust  levels  over  the  existing  Redesigned  Solid  Rocket  Motor  (RSRM) 
and  thus  augment  the  lifting  capacity  of  the  space  shuttle  orbiter.  To  achieve  these 
higher  thrust  levels  and  improve  motor  reliability,  advanced  motor  design  concepts 
were  employed.  In  the  head  end  of  the  motor,  for  instance,  the  propellent  cast  has 
been  changed  from  the  conventional  annular  configuration  to  a "multi-slot" 
configuration  in  order  to  increase  the  burn  surface  area  and  guarantee  rapid  motor 
ignition.  In  addition,  the  igniter  itself  has  been  redesigned  and  currently  features  12 
exhaust  ports  in  order  to  channel  hot  igniter  combustion  gases  into  the 
circumferential  propellent  slots.  Due  to  the  close  proximity  of  the  igniter  ports  to 
the  propellent  surfaces,  new  concerns  over  possible  propellent  deformation  and 
erosive  burning  have  arisen.  The  following  documents  the  effort  undertaken  using 
computational  fluid  dynamics  to  perform  a flow  field  analysis  in  the  top  end  of  the 
ASRM  motor  to  determine  flow  field  properties  necessary  to  permit  a subsequent 
propellent  fin  deformation  analysis  due  to  pressure  loading  and  an  assessment  of 
the  extent  of  erosive  burning. 
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Near  Field  Velocity  Vectors 
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IGNITION  TRANSIENT  CALCULATIONS  IN  THE 
SPACE  SHUTTLE  SOLID  ROCKET  MOTOR 

Rhonald  M.  Jenkins  and  Winfred  A.  Foster,  Jr. 

Aerospace  Engineering  Department 
Auburn  University,  Alabama 


ABSTRACT 

The  work  presented  is  part  of  an  effort  to  develop  a multi- 
dimensional ignition  transient  model  for  large  solid  propellant 
rocket  motors.  On  the  Space  Shuttle,  the  ignition  transient  in  the 
slot  is  induced  when  the  igniter,  itself  a small  rocket  motor,  is 
fired  into  the  head-end  portion  of  the  main  rocket  motor.  The 
computational  results  presented  in  this  paper  consider  two 
different  igniter  configurations.  The  first  configuration  is  a 
simulated  Space  Shuttle  RSRM  igniter  which  has  one  central  nozzle 
that  is  parallel  to  the  centerline  of  the  motor.  The  second 
igniter  configuration  has  a nozzle  which  is  canted  at  an  angle  of 
45°  from  the  centerline  of  the  motor.  This  paper  presents  a 
computational  fluid  dynamic  (CFD)  analyses  of  certain  flow  field 
characteristics  inside  the  solid  propellant  star  grain  slot  of  the 
Space  Shuttle  during  the  ignition  transient  period  of  operation  for 
each  igniter  configuration.  The  majority  of  studies  made  to  date 
regarding  ignition  transient  performance  in  solid  rocket  motors 
have  concluded  that  the  key  parameter  to  be  determined  is  the  heat 
transfer  rate  to  the  propellant  surface  and  hence  the  heat  transfer 
coefficient  between  the  gas  and  the  propellant.  In  this  paper  the 
heat  transfer  coefficients,  pressure  and  velocity  distributions 
are  calculated  in  the  star  slot.  In  order  to  validate  the 
computational  method  and  to  attempt  to  establish  a correlation 
between  the  flow  field  characteristics  and  the  heat  transfer  rates 
a series  of  cold  flow  experimental  investigations  were  conducted. 
The  results  of  these  experiments  show  excellent  qualitative  and 
quantitative  agreement  with  the  pressure  and  velocity  distributions 
obtained  from  the  CFD  analysis.  The  CFD  analysis  utilized  a 
classical  pipe  flow  type  correlation  for  the  heat  transfer  rates. 
The  experimental  results  provide  an  excellent  qualitative 
comparison  with  regard  to  spatial  distribution  of  the  heat  transfer 
rates  as  a function  of  nozzle  configuration  and  igniter  pressure. 
The  results  indicate  that  from  a quantitative  point  of  view  that 
the  pipe  flow  correlation  gives  reasonably  good  results. 
Furthermore,  there  appears  to  be  a direct  correlation  between 
igniter  pressure  and  an  average  Reynolds  number  in  the  star  grain 
slot.  This  may  lead  to  a simple  method  for  modifying  the 
convection  heat  transfer  correlation.  Calculated  results  of 
pressure-vs-time  for  the  first  200  msec  of  motor  firing  of  the 
Space  Shuttle  RSRM  support  the  trends  shown  for  the  heat  transfer 
rate  comparisons  between  the  cold  flow  CFD  and  experimental  data. 
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OUTLINE 


Presentation  of  Cold-Flow  Heat  Transfer/  Aerodynamic 
Results  for  Flow  in  the  Head-End  Star  Slot  of  the  Space 
Shuttle  SRM 

• Single  Port  Igniter 

• 45°  Igniter  (no  center  port) 


Presentation  of  Flow  Visualization  Results  for  Three 
Igniter  Configurations 

• Single  Port  Igniter 

• 45°  Igniter  (no  center  port) 

• 45°  Igniter  (with  center  port) 
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HEAT  TRANSFER  MODEL 


Heat  transfer  from  gas  to  wall  is  assumed  to  be  by 
convection  only. 


The  convection  correlation  utilized  was  originally 
derived  for  turbulent  pipe  flow,  but  can  be  adapted  to 
arbitrary  geometries.  The  correlation  utilized  is  of 
the  form 

_ 0. 152i?e0 -9Pr 

xVU  - — : 

0 . 833  [2 . 2 5 In  (0 . 114Pe0,9)  + 13.2Pr  - 5.8] 

where  Re  is  based  on  the  hydraulic  diameter  of  the 
slot . 
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Calculated  Heat  Transfer,  Single  Port  Igniter 
PlflnlUf  a 1 00  pal 


Measured  Heat  Transfer,  Single  Port  Igniter 
Plgnlttr  s 1 00  psl 
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Calculated  Heat  Transfer,  Single  Port  Igniter 
Plgnlter  a 500  ptl 


Measured  Heat  Transfer,  Single  Port  Igniter 
Pinter  = 500  pal 
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Calculated  Heat  Transfer,  Single  Port  Igniter 
Plgnltor  » 1 500  psl 


Measured  Heat  Transfer,  Single  Port  Igniter 
Pi«nrt*r  *1500  psi 
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Calculated  Heat  Tranafer,  45°  Igniter 
Plgnlter  a 1 00  pal 


Measured  Heat  Transfer,  45°  Igniter 
Plgnitor  *100  psl 
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Calculated  Heat  Transfer,  45°  Igniter 
Plgnlter  = 500  psl 


Measured  Heat  Transfer,  45°  Igniter 
PlgnHif  = 500  psl 
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Calculated  Heat  Transfer,  45°  Igniter 
s 1 500  psl 


Measured  Heat  Transfer,  45°  Igniter 
Piflnit*  = 1 500  psl 
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Igniter  Pressure,  psl 


450  Igniter  (no  center  port)  average  slot  heat  transfer 
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Sensitivity  of  Initial  Shuttle  Head-End  Pressure 
to  Assumed  Heat  Transfer  Correlation 
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Conclusions:  Work-to-Date 


Correlation  between  the  convection  heat  transfer 
model  utilized  and  the  measured  (cold  flow)  values 
of  heat  transfer  is  generally  good. 


For  the  single  port  igniter,  the  heat  transfer 
model  utilized  increasingly  under-predicts  the  heat 
transfer  as  igniter  pressure  increases. 


For  the  45°  igniter  with  no  center  port,  the  heat 
transfer  model  under-predicts  the  heat  transfer 
for  all  igniter  pressures,  with  increasing  under- 
prediction as  igniter  pressure  increases. 


There  appears  to  be  a direct  correlation  between 
igniter  pressure  and  an  average  Reynolds  number 
in  the  star  grain  slot.  This  may  lead  to  a simple 
method  for  modifying  the  convection  heat  transfer 
correlation. 


Calculated  results  of  pressure-vs-time  for  the  first 
200  msec  of  motor  firing  of  the  Space  Shuttle  SRM 
tend  to  support  the  idea  that  heat  transfer  is 
under-predicted  for  higher  igniter  pressures. 
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Abstract 


Under  a NASA  MSFC  SBIR  Phase  I effort  an  interactive  software  package  has  been  developed 
for  the  analysis  of  discrete  (particulate)  phase  dynamics  in  two-phase  flows  in  which  the  discrete 
phase  does  not  significantly  affect  the  continuous  phase.  This  package  contains  a Graphical  User 
Interface  (based  on  the  X Window  system  and  the  Motif" tool  kit)  coupled  to  a particle  tracing 
program,  which  allows  the  user  to  interactively  set  up  and  run  a case  for  which  a continuous  phase 
grid  and  flow  field  are  available.  The  software  has  been  applied  to  a solid  rocket  motor  problem,  to 
demonstrate  its  ease  of  use  and  its  suitability  for  problems  of  engineering  interest,  and  has  been 
delivered  to  NASA  Marshall  Space  Flight  Center. 


I.  Introduction 

Multiphase  flow  effects  in  liquid  and  solid  propulsion  systems  have  profound  implications  on 
performance  and  durability.  For  example,  the  dynamics  of  aluminum  oxide  particulates  in  an  SRM 
affects  the  slag  accumulation,  thereby  affecting  the  performance.  A second  example  is  particulate 
impingement  on  the  motor  casing  and  nozzle  affecting  durability  via  its  influence  on  the  thermal  load 
on  the  insulator.  The  effect  of  the  discrete  phase  on  performance  and  durability  can  be  very 
significant  even  when  the  concentration  of  the  discrete  phase  is  low  enough  to  not  alter  the 
continuous  phase  flow  field  in  a significant  manner.  Under  the  present  effort  a workstation-based 
analysis  tool  has  been  developed  which  can  be  used  by  analysts  and  designers  to  interactively 
assess  the  discrete  phase  effects  during  the  development  and  testing  of  rocket  propulsion  systems. 
The  workstation  based  software  uses  a Lagrangian  analysis  for  the  discrete  phase  motion  and 
includes  a Graphical  User  Interface  allowing  the  user  to  interactively  change  the  relevant  parameters 
to  conduct  parametric  studies.  To  demonstrate  the  suitability  of  the  interactive  software  and  its  ease 
of  use,  it  has  been  successfully  applied  to  a two-dimensional  solid  rocket  motor  test  problem. 

The  software  consists  basically  of  two  parts:  (i)  the  Graphical  User  Interface  for  input/output, 
and  (ii)  the  computational  analysis  for  the  discrete  phase  dynamics.  As  description  of  these  is  given 
below. 
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The  present  analysis  is  based  on  the  Lagrangian  analysis  described  by  de  Jong  et  al.  [1]. 
Although  originally  part  of  a fully  coupled  two-phase  flow  analysis  (see,  for  example,  Madabhushi 
et  al.  [2]  or  de  Jong  and  Sabnis  [3]),  this  analysis  has  been  adapted  for  stand-alone  use,  and 
does  ngt,  in  its  current  form,  depend  on  any  flow  solver.  For  a given  computational  grid  and 

continuous  phase  flow  field,  the  stand-alone  particle  code  integrates  a Lagrangian  equation  of  motion 
for  each  particle.  A brief  description  of  its  relevant  features  is  presented  below. 


The  equation  of  motion  for  a particle  can  be  written  in  the  form 

= r o> 

dt*  m 

where  X is  the  position  vector  of  the  particle,  m is  the  particle  mass,  and  F is  the  total  force  acting  on 
the  particle.  In  general,  F contains  'body'  forces,  such  as  those  due  to  gravity  or  electro-magnetic 
fields,  and  the  force  acted  upon  the  particle  by  the  fluid  (F p).  Integration  of  Eq.  (1)  yields 


-t 


P _ dX 

- dr  + — 
m dt 


At  this  stage  a coordinate  transformation  Y = Y(X)  is  introduced  to  transform  the  equation  of 
motion  into  the  computational  coordinate  space  corresponding  to  the  grid  used  for  the  continuous 
phase  analysis,  i.e. 


yj  = yj(x1/x2/x3) 


(j  = 1,2,3) 


where  yi  are  the  computational  coordinates  (the  components  of  Y)  used  in  the  continuous  phase 
analysis  and  x;  are  the  physical  coordinates  (the  components  of  X).  Let  J be  the  inverse  of  the 
transformation  matrix,  i.e.  let  J be  the  matrix  with  elements  Jjj,  where 

T _ aVL 
JiJ  “ axj 


Then 


dY  dX 

dt  dt 


Substituting  Eq.  (2)  into  Eq.  (3)  yields 


J “ 
dt 
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This  equation  can  be  integrated  (assuming  that  J,  F,  and  m are  constant  over  the  integration  interval 
At)  to  yield 

AY  = %At2J  J + At JUp | (5) 

m v t0 

where  aY  is  the  change  in  the  coordinate  space  position  vector  (i.e.  Y)  during  the  period  At.  Equation 
(5)  then  provides  the  change  in  particle  location  (in  computational  space)  for  a given  time  interval  At 
and  starting  location  corresponding  to  t = \q. 

The  transformation  of  the  Lagrangian  equations  of  motion  from  the  physical  coordinates  to  the 
computational,  body-fitted  coordinates  used  for  the  continuous  phase  analysis  offers  some  very 
significant  advantages.  Since  the  motion  of  a particle  is  tracked  in  the  computational  coordinate 
space  (via  Eq.  (5)),  no  searching  is  needed  for  the  mesh  cell  in  which  a computational  particle 
resides.  This  facilitates  the  computation  of  the  force  on  the  particle  in  the  Lagrangian  equation  of 
motion.  Also,  tracking  the  motion  in  the  computational  space  allows  the  estimation  of  the  time  taken 
by  the  particle  to  cross  the  boundaries  of  the  grid  cell  where  it  currently  resides.  This  estimate  is  used 
in  the  selection  of  sub-time  steps  (described  below)  used  in  integrating  the  equation  of  motion. 
Finally,  the  search  for  boundaries  is  simplified,  which  makes  it  easier  to  determine  whether  a particle 
has  reached  a solid  boundary  and  should  be  reflected,  or  a particle  has  left  the  computational 
domain  and  should  be  omitted.  This  latter  advantage  is  of  major  importance  when  the  number  of 
boundaries  is  large  (see  de  Jong  et  al.  [1]  and  McConnaughey  et  al.  [4]  for  an  interesting 
application). 

Each  particle  is  moved  through  the  domain  by  integrating  its  Lagrangian  equation  of  motion 
using  a sequence  of  "sub-time  steps'  Ats.  These  sub-time  steps  are  chosen  as  the  minimum  time 
for  the  particle  to  cross  the  nearest  cell  boundary  and  a kinematic  time  scale  defined  as  a fraction 
(e.g.  0.1)  of  the  ratio  of  particle  velocity  to  acceleration.  It  should  be  noted  that,  for  accuracy 
reasons,  the  maximum  allowable  sub-time  step  could  be  significantly  different  for  different 
particles  depending  upon  their  location,  mass,  etc.  The  sub-time  step  approach  has  been  adopted 
to  provide  an  accurate  representation  of  the  particle  motion,  since  at  the  end  of  each  sub-time 
step  the  force  on  the  particle  is  updated  using  the  particle  attributes  and  the  continuous  phase 
values  corresponding  to  the  new  particle  locations. 


11.2  Force  on  a Particle 


In  Equation  (1),  the  force  F on  the  particle  has  been  assumed  to  consist  of  a drag  force,  a 
gravitational  acceleration  force,  and  a force  due  to  the  pressure  gradient  in  the  liquid, 

P = Fn  + mg  — — yp  (6) 

p P 

where  g is  the  (gravitational)  acceleration  vector,  Pp  is  the  particle  density,  p is  the  local  pressure  in 
the  continuous  phase  and  where  the  drag  force  Fq  is  given  by 

| cD'”'Dp2lDRl0R  <7> 

Here  p is  the  density  of  the  continuous  phase,  Dp  is  the  particle  diameter  and  Ur  is  the  velocity  of 
the  particle  relative  to  the  continuous  phase  velocity.  Cq  is  the  drag  coefficient,  which  is  related 
to  the  particle  Reynolds  number 
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Re„  = 


pURDp 


'P  n 

(where  fi  is  the  viscosity  of  the  continuous  phase)  via  the  expression 


CD  - 


24  fl  + 1 Re„2/3] 

for 

Ren 

< 1000 

Rep  l 6 P J 

P 

0.424 

for 

Rep 

> 1000 

(9) 


The  second  term  on  the  right  hand  side  of  Eq.  (6)  represents  the  gravitational  (or  other  acceleration) 
effects,  which  could  be  significant  for  solid  particles  or  liquid  droplets  in  a gaseous  continuous  phase. 
The  effect  of  buoyancy,  on  the  other  hand,  which  is  included  in  the  last  term  of  Eq.  (6),  may  not 
always  be  of  importance. 


11.3  Discrete  Phase  Turbulent  Dispersion 

In  two-phase  flows,  the  continuous  phase  turbulence  results  in  dispersion  of  the  discrete  phase. 
This  process  is  modeled  by  some  researchers  (see  Abbas  et  al.  [5])  by  adding  a "diffusion"  velocity 
obtained  from  some  phenomenological  model  to  the  particle  velocity.  In  the  present  analysis,  this 
effect  is  modeled  by  evaluating  the  fluid  force  on  a particle  using  an  instantaneous  continuous  phase 
velocity  field  rather  than  a mean  velocity  field.  The  instantaneous  velocity  components  are  obtained 
by  adding  stochastically  generated  turbulent  velocity  components  to  the  mean  velocity  field  for  the 
continuous  phase.  If  the  continuous  phase  turbulence  is  assumed  to  be  isotropic,  and  the  random 
turbulent  velocity  components  are  assumed  to  have  a Gaussian  probability  distribution  with  standard 
deviation  a,  then 

- - ( ? r 


where  k is  the  turbulent  kinetic  energy,  and  the  random  turbulent  velocity  components  can  be 
obtained  as 


4k 

3 


) 


1/2 


erf _1 (2x  - 1) 


(10) 


where  x is  random  variable  with  uniform  probability  distribution  between  0 and  1 (i.e.  x is  a random 
number  between  0 and  1).  This  latter  result  follows  from  the  observation  that  if  x is  random  variable 
with  uniform  probability  distribution  between  0 and  1,  then  F'1  (x)  is  a random  variable  whose 
cumulative  probability  function  is  F.  This  technique  offers  a way  to  treat  the  effect  of  gas  phase 
turbulence  on  the  particulate  phase  dispersion  in  a manner  which  is  less  empirical  than  the 
techniques  traditionally  used  with  Eulerian-Eulerian  analyses.  Similar  techniques  for  modeling  the 
discrete  phase  turbulent  dispersion  have  been  used  by,  for  example,  Dukowicz  [6],  Gosman  and 
loannides  [7],  and  Hotchkiss  [8]. 

If  the  turbulent  energy  k is  not  directly  available,  but  an  eddy  viscosity  nj  and  a mixing  length  A 
are  (for  example,  when  a mixing  length  turbulence  model  is  used  in  the  continuous  phase  analysis), 
then  k can  be  obtained  from  the  relation 


\ 

pk2f 


(12) 
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where  is  a function  of  p j that  can  be  determined  from  the  expression 

2.5 


= 0.09  exp 


Mt 

1 + 0.02  — C,, 


(13) 


(cf.  Launder  and  Spalding  [10]). 


III.  Graphical  User  Interface 

The  Graphical  User  Interface  (GUI)  developed  under  the  present  effort  is  based  on  the  X 
Window  system  and  the  Motif  tool  kit.  Some  advantages  to  using  X and  Motif  include  improved 
code  portability,  maintaining  a consistent  look  and  feel  between  applications  and  the  ability  to  run  the 
interface  on  one  machine  and  to  display  it  on  another.  Its  portability  was  verified  by  compiling  and 
running  the  GUI  both  on  SGI  4D/25TG  and  IBM  RS/6000  Workstations.  The  main  window  of  the  GUI 
contains  a menubar  with  several  pull-down  menus  (Fig.  1)  that  allow  the  user  to  perform  functions 
associated  with  file  I/O,  problem  specification,  or  output  display.  In  most  cases,  selection  of  a menu 
item  from  a pull-down  menu  opens  up  a property  window  that  provides  the  user  with  the  requested 
input  items.  Depending  on  the  item,  user  input  is  possible  via  type-ins,  toggle  buttons,  selection  from 
lists,  or  graphically  by  using  the  mouse.  Some  of  the  features  of  the  GUI  are: 

(i)  Grid  and  flow  field  files  can  be  used  in  standard  PLOT3D-type  format. 

(ii)  Boundaries  can  be  defined  in  the  computational  domain  by  adding  or  deleting  ■zones'  or 
■blocks'.  Here  a "zone’  is  defined  as  a region  inside  which  there  is  a flow  field,  while  a 'block'  is 
defined  as  a solid  region  (without  a flow  field).  These  "zones'  or  'blocks'  can  be  specified 
graphically  by  using  the  mouse  (Figs.  2-3).  This  approach  allows  the  user  to  define  complex 
(2-D)  geometries  without  the  need  for  code  modifications. 

(iii)  To  define  the  boundary  type  (such  as  inlet/exit,  wall,  or  symmetry  line)  and  the  properties 
(restitution  and  friction  coefficients)  of  a solid  wall,  this  boundary  can  be  selected  by  "clicking" 
on  it.  The  boundary  type  can  then  be  set  by  choosing  the  appropriate  one  from  a menu  (Fig. 
4),  while  the  relevant  properties  can  be  typed  in. 

(iv)  The  injection  locations  can  be  specified  graphically  by  clicking  on  the  desired  grid  points  in 
the  computational  domain  (Fig.  5).  The  particle  properties  at  these  locations  (such  as  the 
particle  velocity,  its  radius  and  its  density)  can  be  specified  via  type-ins.  By  defining  default 
properties,  only  those  properties  that  change  from  one  particle  to  the  next  have  to  be  typed  in. 

(v)  Wherever  possible,  error  and  consistency  checks  are  performed  by  the  interface  on  user  input 
prior  to  running  the  particle  code.  This  helps  to  produce  a more  accurate  input  file  and 
prevents  needless  execution  of  the  particle  code. 

The  Graphical  User  Interface  is  coupled  to  the  particle  code,  so  that  it  can  be  used  both  to  set 
up  and  to  run  a particular  problem.  Execution  of  the  particle  tracking  program  under  the  GUI  displays 
the  particle  traces  while  they  are  being  computed;  once  the  program  is  finished,  the  image  can  be 
manipulated  interactively.  A complete  description  of  all  the  functions  included  in  the  GUI  and 
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guidelines  for  setting  up  and  running  a new  case  can  be  found  in  the  User’s  Manual;  on-line  help  is 
available  for  most  functions. 


IV.  Application  and  Results 

To  demonstrate  the  suitability  of  the  software  package  developed  under  the  present  effort  and 
its  ease  of  use  for  solid  rocket  motor  applications,  the  Particle  Trace  Interface  was  used  to  set  up  and 
run  particle  trajectories  in  the  so-called  super  BATES  motor.  The  geometry  of  this  motor  is  shown  in 
Fig.  6.  The  flow  field  in  this  motor  was  computed  with  SRA’s  MINT  code,  on  the  grid  shown  in  Fig.  6. 
This  grid  contains  two  ■cut-out"  regions,  adjacent  to  the  grain  surface  and  the  nozzle  wall.  The  "full* 
grid  is  shown  in  Fig.  7.  This  grid,  and  the  flow  field,  were  then  used  as  input  to  the  Particle  Trace 
Interface  (PTI).  Using  the  PTI,  the  boundaries  of  the  computational  domain  were  constructed  by 
putting  two  "blocks’  into  the  "zone"  that  comprised  the  full  grid  (cf.  Fig.  5).  Alternatively,  it  would  have 
been  possible  to  build  up  the  domain  as  the  union  of  three  "zones",  corresponding  to  the  combustion 
chamber,  the  slot  region,  and  the  nozzle  region.  At  this  point  it  should  be  noted  that  Fig.  7 was 
actually  obtained  by  displaying  the  grid  in  the  PTI  after  it  had  been  read  in,  while  Fig.  6 was  obtained 
by  displaying  the  grid  after  the  boundaries  had  been  specified.  Next,  the  PTI  was  used  to  define  the 
"boundary  conditions'  (insofar  as  these  are  needed  to  determine  the  behavior  of  a particle  when  it 
reaches  a boundary)  and  a set  of  particle  injection  locations  (the  solid  circles  in  Fig.  5).  Particle  and 
flow  field  information  was  specified  to  complete  the  problem  setup,  after  which  the  PTI  was  used  to 
(interactively)  execute  the  particle  tracing  code.  Figs.  8 and  9 show  the  resulting  particle  traces 
(again  obtained  from  the  PTI).  The  whole  operation  was  completed  without  exiting  from  the  PTI. 
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Figure  1.  Main  Panel  of  the  Particle  Trace  Interface  with  the  File  Pulldown  Menu. 
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Figure  2.  The  Boundary  Definition  Window  with  the  Blocks  Pulldown  Menu. 
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Figure  3.  Adding  a Block. 
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Figure  4.  The  Boundary  Type  Selection  Window. 
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Figure  6.  Super  BATES  Motor  Geometry  an  Grid  in  the  Physical  Domain 
with  Blocked-out  Regions. 
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Figure  7.  Full  Grid  in  the  Physical  Domain. 
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Figures.  Particle  Trajectories. 
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Figure  9.  Particle  Traces  Overlaid  on  the  Grid; 

Enlarged  View  near  the  Aft  End  of  the  Motor. 
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Statement  of  Problem 

There  is  a critical  shortage  of  data  pertaining  to  the  flowfield  characteristics  in  a liquid  propellant 
rocket  chamber  for  hot-fire  conditions.  For  a liquid  oxygen  (LOX)/gaseous  hydrogen  (GHj)  propellant 
combination,  either  shear  or  swirl  coaxial  injectors  are  typically  used  to  atomize  the  liquid  propellant  into 
drops.  Understanding  the  atomization  process  under  hot-fire  conditions  represents  the  first  step  in 
understanding  the  subsequent  processes  of  vaporization,  mixing  and  combustion.  The  flowfield  here  is 
two-phase  and  therefore  experiments  that  detail  the  intact  liquid  jet,  drop  size  and  velocity,  and 
combustion  length  are  necessary  for  understanding  the  physics  of  the  problem. 

Objective  of  Work 

The  objective  of  the  current  work  is  to  experimentally  characterize  the  flowfield  associated  with 
an  uni-element  shear  coaxial  injector  burning  LOX/GH2  propellants.  These  experiments  were  carried  out 
in  an  optically-accessible  rocket  chamber  operating  at  a high  pressure  (*400  psia).  Quantitative 
measurements  of  drop  size  and  velocity  were  obtained  along  with  qualitative  measurements  of  the 
disintegrating  jet. 

i 

Approach 

The  experiments  were  conducted  at  the  Cryogenic  Combustion  Laboratory  at  Penn  State 
University.  This  laboratory  provides  the  capability  for  firing  both  gaseous  and  liquid  propellant  sub-scale 
rocket  engines.  A modular  rocket  chamber  which  provides  extensive  optical  access  was  designed  for  the 
experiments.  The  cross-section  of  the  rocket  is  50.8  mm  (2  in.)  square  and  the  chamber  length  which 
can  be  easily  varied  was  248  mm  (9.75  in.).  The  flowfield  downstream  of  a shear  coaxial  injector  was 
characterized  using  laser-based  diagnostic  techniques.  The  inner  diameter  of  the  LOX  post  was  3.43  mm 
(0.135  in.)  and  the  post  was  recessed  3.78  mm  (0.15  in.).  The  inner  diameter  of  the  fuel  annulus  was 
4.19  mm  (0.165  in.)  and  the  outer  diameter  was  7.11  mm  (0.28  in.).  The  nominal  LOX  and  GH2 
flowrates  were  0.1 1 kg/s  (0.25  Ibm/s)  and  0.021  kg/s  (0.047  lbm/s)  respectively,  thus  resulting  in  a 
nominal  O/F  ratio  of  5.3:1.  These  flow  rates,  coupled  with  the  nozzle  dimensions  yielded  a chamber 
pressure  of  2.74  MPa  (*400  psia).  The  duration  of  a test  run  was  four  seconds. 

A Phase  Doppler  Particle  Analyzer  (PDPA)  was  used  to  measure  LOX  drop  size  and  velocity  at 
various  radial  locations  for  an  axial  position  63.5  mm  (2.5  in.)  downstream  of  the  injector  face.  The 
measured  Sauter  Mean  diameter  (SMD)  ranged  from  1 10  jun  at  the  centerline  to  about  60  /zm  at  a 
9.5  mm  (0.375  in.)  radial  location.  At  greater  radial  locations,  no  drops  were  observed.  The  results 
indicate  that  under  hot-fire  conditions,  the  drops  formed  from  the  shear  coaxial  injector  are  confined  to 
a narrow  circumferential  region. 

Conclusions 

These  experiments  represent  the  first  time  that  drop  sizes  have  been  measured  under  combusting 
conditions  for  cryogenic  propellants.  These  results  are,  in  general,  encouraging  with  respect  to 
applications  of  laser-based  diagnostics  to  LOX/GH-,  uni-element  rocket  studies.  A comprehensive 
mapping  of  the  flowfield  will  need  to  be  completed  to  gain  a thorough  understanding  of  the  physics  of 
this  complex  problem. 
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MOTIVATION 


To  obtain  fundamental  data  under  well  characterized 
conditions 

- For  gaseous  propellant  rockets 

- For  liquid  propellant  rockets 

- For  various  injector  types 

Fundamental  data  would  form  the  basis  for 

- Empirical  correlation  validation 

- CFD  code  validation 
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BACKGROUND 


Atomization  models  typically  are: 

- Anchored  to  cold  flow  experimental  results 

We  and  Re  of  cold  flow  experiments  differ  by  an 
order  of  magnitude  from  actual  rocket  conditions 

Results  have  to  be  extrapolated 

- Predicted  from  analytical  models  based  on  linear 
stability  theory 


Drop  size  data  for  hot-fire  conditions  would: 

- Validate  atomization  models 

- Validate  methodology  for  extending  cold  flow  data  to 
hot  flow  conditions 

- Develop  hot-flow  correlations  for  direct  use  in 
combustion  models 
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OBJECTIVE 


To  characterize  the  flowfield  downstream  of  a shear 
coaxial  injector  using  LOX/GH2  propellants  under 
combusting  conditions 

- Drop  size  and  velocity  measurements  using  Phase 
Doppler  Interferometry 

- Laser  sheet  visualizations  of  near  breakup  region 
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FACILITY  CAPABILITIES 


Propulsion  Engineering  Research  Center 
Cryogenic  Combustion  Laboratory 


Propellants: 

Gaseous  Hydrogen 
Gaseous  Methane 
Gaseous  Oxygen 
Liquid  Oxygen 

Flow  rates  (maximum): 

Gaseous  oxygen:  0.1  lbm/s 

Liquid  oxygen:  1.0  lbm/s 

Gaseous  hydrogen:  0.25  lbm/s 


Typical  mixture  ratios:  4-8 

Maximum  chamber  operating  pressure:  1000  psi 

Typical  test  time:  1 - 5 s 

Longer  tests  subject  to  hardware  cooling  and  gas 
supply  specifications 
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PENN  STATE  ROCKET 
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Cross-Sectional  View  of  the  Windowed  Rocket  Chamber 
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TEST  CONDITIONS 


Chamber  pressure 
LOX  flowrate 
GH2  flowrate 
Mixture  ratio  (O/F) 
LOX  velocity 
GH2  velocity 
Momentum  ratio  (F/O) 
Velocity  ratio  (F/O) 

Re1 

We2 

lRe=p,U,  d/p, 

2 W e = pg(Ug-U,)2d/  a 


2.67  Mpa  (387  psia) 
0.112  kg/s  (0.245  lbm/s) 
0.021  kg/s  (0.045  lbm/s) 
5.4 

13.5  af/s  (44.1  ft/s) 

381  m/s  (1250  ft/s) 

5.22 

28.3 

5.03x10s 

2.06xl05 
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ROCKET  PARAMETERS 
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PDPA  RESULTS 
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SHEAR-COAXIAL  INJECTOR  SPRAYS 
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Weg  (=pg(Ug-U1)2d/ff) 


DROP  SIZE  EQUATIONS 


Ui 


< 

< 

CO 

< 

z 

1 

2 
< 
> 

2 


Q. 

•5 


rls- 

a 


& 


isf 


Uj  co 


ORIGIN  a;.  ^ 

Of  POOR  QUALITY 


631 


DROP  SIZE  PREDICTIONS 
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HOT  FIRE/COLD  FLOW  COMPARISON 
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(Ug-U,)2d/a)  2.06  x 10 


HOT  FIRE/COLD  FLOW  COMPARISON 

Sauter  Mean  Diameter  (D32) 
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HOT  FIRE/COLD  FLOW  COMPARISON 
Mean  Drop  Velocity  (UD) 
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SUMMARY 


• Measured  liquid  oxygen  drop  size  and  velocity  in 
combusting  environment 

- Intact  core  extends  well  beyond  the  injector 

- Drops  confined  to  narrow  region 

• Correlations  based  on  cold  flow  data  inadequate  for 
predicting  drop  size  in  LOX/GH2  combusting  flow 

• Compared  drop  measurements  between  cold  flow  and 
combusting  conditions  for  similar  liquid  and  gas 
flowrates 

- Re,  and  Weg  differ  by  an  order  of  magnitude 

- Mean  drop  size  larger  for  hot  fire  conditions 
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by 

Edward  J.  Mularz 

U.  S.  Army  Vehicle  Propulsion  Directorate,  ARL 
NASA  Lewis  Research  Center 

Daniel  L.  Bulzan 
NASA  Lewis  Research  Center 

Kuo-Huey  Chen 
The  University  of  Toledo 


ABSTRACT 

The  next  generation  of  commercial  aircraft  will  include  turbofan  engines  with 
performances  significantly  better  than  those  in  the  current  fleet.  Control  of  particulate  and 
gaseous  emissions  will  also  be  an  integral  part  of  the  engine  design  criteria.  These  performance 
and  emission  requirements  present  a technical  challenge  for  the  combustor:  control  of  the  fuel 
and  air  mixing  and  control  of  the  local  stoichiometry  will  have  to  be  maintained  much  more 
rigorously  than  combustors  in  current  production.  A better  understanding  of  the  flow  physics  of 
liquid  fuel  spray  combustion  is  necessary.  This  presentation  describes  recent  experiments  on 
spray  combustion  where  detailed  measurements  of  the  spray  characteristics  were  made,  including 
local  drop-size  distributions  and  velocities.  In  addition,  an  advanced  combustor  CFD  code  has 
been  under  development  and  predictions  from  this  code  are  presented  and  compared  with 
measurements.  Studies  such  as  these  will  provide  information  to  the  advanced  combustor 
designer  on  fuel  spray  quality  and  mixing  effectiveness.  Validation  of  new  fast,  robust,  and 
efficient  CFD  codes  will  also  enable  the  combustor  designer  to  use  them  as  valuable  design  tools 
for  optimization  of  combustor  concepts  for  the  next  generation  of  aircraft  engines. 
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COMBUSTOR 
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Coflow 


Gas  Phase  Velocity  at  5 mm  Downstream 

Mean  Axial  Velocity 
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Gas  Phase  Velocity  at  5 mm  Downstream 

Fluctuating  Axial  Velocity 
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NUMERICAL  ALGORITHM 
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Difficulties  with  Compressible  Flow  Algorithms  at  Low  Mach 
Numbers 
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METHOD  OF  APPROACH 
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Pressure  rescaled  so  that  all  eigenvalues  have  the  same  order  of  mag- 
nitude. Physical  acoustic  waves  removed  and  replaced  with  pseudo- 
acoustic waves  which  travel  at  speed  comparable  to  fluid  convective 


Mean  Velocity  Profiles 
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Velocity  Vectors 
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Preliminary  Results  Obtained  for  Two-Phase, 
Turbulent,  Reacting  Combustor  Flowfield 
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PROGRESS  IN  ADVANCED  SPRAY  COMBUSTION 
CODE  INTEGRATION 

Pak-Yan  Liang 

Rocketdyne  Division,  Rockwell  International 
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ABSTRACT 

A multiyear  project  to  assemble  a robust,  multiphase  spray  combustion  code  is 
now  underway  and  gradually  building  up  to  full  speed.  The  overall  effort 
involves  several  university  and  government  research  teams  as  well  as 
Rocketdyne.  The  first  part  of  this  paper  will  give  an  overview  of  the  respective 
roles  of  the  different  participants  involved,  the  master  strategy,  the  evolutionary 
milestones,  and  an  assessment  of  the  state-of-the-art  of  various  key 
components.  The  second  half  of  this  paper  will  highlight  the  progress  made  to- 
date  in  extending  the  baseline  Navier-Stokes  solver  to  handle  multiphase, 
multispecies,  chemically  reactive  sub-  to  supersonic  flows.  The  major  hurdles 
to  overcome  in  order  to  achieve  significant  speed  ups  are  delineated  and  the 
approaches  to  overcoming  them  will  be  discussed. 
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SEVERAL  PHYSICAL  PROCESSES  IN  LIQUID 
PROPULSION:  I.E.,  ATOMIZATION,  EVAPORATION, 
DENSE  SPRAY  EFFECTS 
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CF.  EXTERNAL  FLOW  AERODYNAMIC  CODES  CIRCA  1980. 
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AS  BASE 

INCORPORATION  OF  NOVEL  ENHANCE- 
MENT FEATURES  CURRENTLY  BEING 
DEVELOPED  IN  OTHER  TECHNOLOGY 
EFFORTS 


APPLICATIONAL  OBJECTIVE 

OPTIMAL  ENGINEER  - MODEL  INTERFACE  (DESIGNER) 
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SYNOPSIS  OF  SOLUTION  APPROACH 
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SOLVE  WITH  STONE’S  STRONGLY  IMPLICIT  PROCEDURE 
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8.  1/2  02  = O (O)  = (KO)  (^02) 

FOR  KINETICS  SPECIES,  PREVIOUS  GENERAL  KINETICS  MODEL 
IS  RETAINED 
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CFD  Analysis  of  Spray  Combustion  and  Radiation 
in  OMV  Thrust  Chamber 


hy 

M.G.  Giridharan,  A.  Krishnan  and  A.J.  Przekwas 
CFD  Research  Corporation 
Huntsville  AL  35805 
and 

K.  Gross 

NASA  Marshall  Space  Flight  Center 
Alabama  35812 


ABSTRACT 

The  Variable  Thrust  Engine  (VTE),  developed  by  TRW,  for  the  Orbit  Maneuvering  Vehicle  (OMV)  uses 
a hypergolic  propellant  combination  of  Monomethyl  Hydrazine  (MMH)  and  Nitrogen  Tetroxide  (NTO) 
as  fuel  and  oxidizer,  respectively.  The  propellants  are  pressure  fed  into  the  combustion  chamber 
through  a single  pintle  infection  element  The  performance  of  this  engine  is  dependent  on  the  pintle 
geometry  and  a number  of  complex  physical  phenomena  and  their  mutual  interactions.  The  most 
important  among  these  are;  (1)  atomization  of  the  liquid  jets  into  fine  droplets;  (2)  the  motion  of  these 
droplets  in  the  gas  held;  (3)  vaporization  of  the  droplets;  (4)  turbulent  mixing  of  the  fuel  and  oxidizer; 
and  (5)  hypergolic  reaction  between  MMH  and  NTO. 

Each  of  the  above  phenomena  by  itself  poses  a considerable  challenge  to  the  technical  community.  In  a 
reactive  flow  field  of  the  kind  occurring  inside  the  VTE,  the  mutual  interactions  between  these 
physical  processes  tend  to  further  complicate  the  analysis. 

The  objective  of  this  work  is  to  develop  a comprehensive  mathematical  modeling  methodology  to 
analyze  the  flow  field  within  the  VTE.  Using  this  model,  the  effect  of  flow  parameters  on  various 
physical  processes  such  as  atomization,  spray  dynamics,  combustion,  and  radiation  is  studied.  This 
information  can  then  be  used  to  optimize  design  parameters  and  thus  improve  the  performance  of  the 
engine. 

The  REFLEQS  CFD  Code  is  used  for  solving  the  fluid  dynamic  equations.  The  spray  dynamics  is 
modeled  using  the  Eulerian-Lagrangian  approach.  The  discrete  ordinate  method  with  12  ordinate 
directions  is  used  to  predict  the  radiative  heat  transfer  in  the  OMV  combustion  chamber,  nozzle,  and 
the  heat  shielcL  The  hypergolic  reaction  between  MMH  and  NTO  is  predicted  using  an  equilibrium 
chemistry  model  with  13  species. 

The  results  indicate  that  mixing  and  combustion  is  very  sensitive  to  the  droplet  size.  Smaller  droplets 
evaporate  faster  than  bigger  droplets,  leading  to  a well  mixed  zone  in  the  combustion  chamber.  The 
radiative  heat  flux  at  combustion  chamber  and  nozzle  wails  are  an  order  of  magnitude  less  than  the 
conductive  heat  flux  Simulations  performed  with  the  heat  shield  show  that  a negligible  amount  of 
fluid  is  entrained  into  the  heat  shield  region.  However,  the  heat  shield  is  shown  to  be  effective  in 
protecting  the  OMV  structure  surrounding  the  engine  from  the  radiated  heat. 
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Estimate  Radiative  Heat  Flux 
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TEMPERATURE  DISTRIBUTION  IN  VTE  CHAMBER/NOZZLE 

(GAS-GAS  MODEL) 
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TEMPERATURE  DISTRIBUTION  IN  VTE  CHAMBER/NOZZLE 

(SPRAY  MODEL) 
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Mach  Number  Distribution  in  VTE  Chamber/Nozzle 

(Spray  Model) 
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EVAPORATED  NTO  MASS  DISTRIBUTION  IN  VTE  CHAMBER 
(100%  POWER  LEVEL) 
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EFFECT  OF  RADIATION  IN  VTE  CHAMBER/NOZZLE 
(TEMPERATURE  DISTRIBUTION,  100%  POWER  LEVEL) 


NO  RADIATION 


WITH  RADIATION 
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(Gas-Gas  Model,  100%  Power  Level) 


EFFECT  OF  RADIATION  IN  VTE  CHAMBER/NOZZLE/HEAT  SHIELD 
(TEMPERATURE  DISTRIBUTION,  100%  POWER  LEVEL) 
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RADIATIVE  FLUX  DISTRIBUTION  IN  VTE  CHAMBER/NOZZLE 

(100%  POWER  LEVEL) 
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Need  an  Advanced  Model  for  Atomization  Due  to 
Impingement  of  Unlike  Liquids 


N95-  23646 


Development  of  an  Atomization  Methodology  for 
Spray  Combustion 


S.  P.  Seung,  C.  P.  Chen 

Department  of  Chemical  Engineering  I 

University  of  Alabama  in  Huntsville 
Huntsville,  AL  35899 

and 

Y.  S.  Chen 

Engineering  Sciences  Inc. 


ABSTRACT 

In  liquid  rocket  propulsion,  the  knowledge  and  the  understanding  of  liquid-gas  interfacial 
phenomena  are  very  important.  This  is  of  keen  importance  for  predicting  the  onset  of 
cavitation  occurring  in  swirl  injection  elements  used  in  STME,  as  well  as  atomization 
processes  in  shear-induced  injector’s  (co-axial)  and  impinging  injector  elements.  From  the 
fact  that  all  the  physical  processes  including  droplets  size  distribution,  droplet  dispersion, 
mixing  and  combustion  are  controlled  by  atomization  processes,  it  is  expected  that  the 
successful  incorporation  of  the  volume  of  fraction  (VOF)  will  greatly  enhance  the  analytical 
capability  for  predicting  spray  combustion  processes  in  liquid-fueled  engines. 

In  this  paper,  a methodology  is  developed  to  define  and  track  interfaces  between  two 
fluids  in  a non-orthogonal, body-fitted  grids  using  a single  fractional  volume  of  fluid(VOF) 
variable  to  describe  the  distribution  of  the  liquid  phase  in  a gas-liquid  flow  field.  This 
method  was  implemented  in  a matured  CFD  code  MAST  (Multiphase  All-Speed  Transient) 
utilizing  the  general  PISO-C  algorithm.  For  the  preliminary  study  for  analyzing  the  spray 
combustion  and  tracking  the  interface  between  two  phase,  we  will  report  the  progress  on 
simulation  of  the  instability  on  the  liquid  column,  the  surface  wave  instability  and  the 
droplet  breakup  from  the  liquid  surface. 
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Coaxial  Atomization  Model  and  Phenomena 
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# SHARP  SURFACE  BETWEEN  TWO  FT  .1 ITDS 
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# SURFACE  TENSION  EFFECT. 
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ree  surface  curvature 


Predictor  Step 
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• First  Corrector  Step 
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# VALIDATIONS  & APPLICATIONS 
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Sauter  mean  diameter  ver.  distance  from  the  injector 
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Test  Conditions  for  the  Measurement  of  Yokota  et  al. 
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o Liquid  Column  Instability  Problem. 
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Liquid  Column  Jet  Breakup  Problem 
Time  = 2.501  x 10‘2  sec  Cycle=  1114 

(a)  Surface 

(b)  Velocity  vector 


Liquid  Column  Jet  Breakup  Problem 
Time  = 3.700  x 10'2  sec  Cycle=  1647 

(a)  Surface 

(b)  Velocity  vector 


Liquid  Column  Jet  Breakup  Problem 
Time  =4.002  x 10’2  sec  Cycle=  1813 

(a)  Surface 

(b)  Velocity  vector 
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743 


Liquid  Droplet  Breakup  Problem 

(a)  Time  = 4.078  x 10'6  sec  Cycle=  26 

(b)  Time  = 1.909  x 10'5  sec  Cycle=  123 


Liquid  Droplet  Breakup  Problem 

(a)  Time  = 3.410  x 10‘5  sec  Cycle=  219 

(b)  Time  = 4.401  x 10'5  sec  Cycle=  284 


Liquid  Droplet  Breakup  Problem 

(a)  Time  = 5.409  x 10'5  sec  Cycle=  351 

(b)  Time  = 5.91 1 x 10  5 sec  Cycle=  385 
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N95- 23647 


MODELING  OF  NON-SPHERICAL  DROPLET  DYNAMICS 

Z.  T.  Deng,  G.  S.  Liaw 

Alabama  A&M  University 
Huntsville,  AL  35762 
(205)  851-5565 

L.  Chou 

Induced  Environment  Branch 
Fluid  Dynamics  Division,  NASA/MSFC 
Huntsville,  AL  35812 


ABSTRACT 

A two-dimensional  time-dependent  computer  code  based  on  the  modified  Arbitrary  La- 
grangian  Eulerian(ALE)  technique  has  been  developed  to  simulate  non-spherical  droplet  dy- 
namics and  evaporation  under  convective  flows  at  real  rocket  combustion  chamber  conditions. 
The  equations  of  mass,  momentum,  energy  and  species  are  simultaneously  solved  for  both 
liquid  and  gas  phases  with  an  accurate  dynamic  interface  tracking.  The  jump  boundary  con- 
ditions across  the  deforming  droplet  surface  are  obtained  by  applying  the  integral  forms  of 
conservation  of  mass,  momentum,  and  energy.  At  each  time  step,  the  interface  geometry  and 
flow  properties  at  the  droplet  surface  are  implicitly  solved  by  satisfying  the  interface  boundary 
conditions.  A Lagrangian  technique  was  developed  to  track  the  arbitrarily  moving  interface 
between  the  liquid  droplet  and  the  external  gas.  An  elliptic  grid  generator  is  adopted  to 
dynamically  reconstruct  grids  both  inside  and  outside  the  droplet  surface. 

This  code  has  been  used  to  study  droplet  oscillation,  droplet  deformation/breakup,  non- 
spherical  droplet  evaporation  in  both  low  and  high  pressure  convective  flows. 

This  presentation  briefly  describes  the  numerical  algorithm  for  modeling  of  the  non- 
spherical  droplet  dynamics  and  demonstrates  the  representative  simulation  results  of  non- 
spherical  droplet  evaporation  at  low  and  high  pressure  convective  flows.  Potential  applications 
of  this  code  to  rocket  combustor  design  and  performance  predictions  are  discussed. 
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Modeling  of  Non-Spherical  Droplet  Dynamics 


I 


^ >. 
d +* 

3 e g 

” s s 

M.5  » 
fl  fl  « 

< 

S3  < S3 

a * 

a 

d 

x 

d 


0) 

Q 


o 

£ 

0) 

XI 

N 


co 

'S 

d 

a 


x\ 

V 

d 
d 
m 

a 

d 'g 
o d 
0) 


XI 

U 

U 


d 

o 


CO 

> 
• H 


CSI 

rH 

X 

IO 

CO 

X 

◄ 

CN 

0) 


P .73 


d 
o 

M 

, «*H 

d > 
d d 
>>& 

* T3 
0) 
O 
d 
'd 
d 


CO 

w 


& 

13 

d 

a 


s 

d 
d 

P ^ 

s ◄ 

M Z 


Pm 


O 

Pm 


750 


1993  Workshop  for  CFD  Application  in  Rocket  Propulsion 


Motivation: 
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in  realistic  combustion  chamber. 


Model  Description: 
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Modified  Arbitrary-Lagrangian-Eulerian  numerical  method. 


Gas-Liquid  Interface  Conditions: 
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s-Liquid  Interface  Conditions  (Cont.): 
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fraction  Ygf  and  surface  pressure  Pg , Pi . 


Gas-Liquid  Interface  Conditions  (Cont.): 
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Radii  of  curvature  at  interface  - Cubic  Spline 


Dynamic  Interface  Tracking: 
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Low-Pressure  Droplet  Evaporation: 
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Parachute  Breakup 
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Gas  Temperature,  Parachute 
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High-Pressure  Droplet  Evaporati 
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HP,  Gas  Velocity,  2M/S 
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Abstract  for  the  Elventh  CFD  Working  Group  Meeting: 

A Fine-Grid  Model  for  the 
ASRM  Aft  Segment  with  Gimballed  Nozzle 

Presented  By:  Dr.  Edward  J.  Reske 

Results  from  computational  fluid  dynamic  analyses  for  complex 
three-dimensional  internal  flows  in  the  Advanced  Solid  Rocket  Motor 
(ASRM)  are  presented.  In  particular,  flow  visualization  and  tabulated 
results  from  a fine-grid  model  consisting  of  1.5  M grid  points  for  the 
ASRM  Aft  Segment  at  the  19-second  burn  time  with  an  8-degree  nozzle 
gimbal  angle  are  shown.  The  results  from  this  model  will  enable 
characterization  of  various  aspects  of  the  ASRM  internal  environment, 
and  in  particular  will  allow  an  assessment  of  the  heat  transfer  and 
stresses  exerted  on  the  submerged  nozzle,  casing  insulation,  and  nozzle- 
case  joint. 
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Structures  and  Dynamics  Laboratory 
Computational  Fluid  Dynamics  Branch 
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Animation  by  Catherine  Dumas, 
Sverdrup  Technology 
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Application  of  CFD  Analyses  to  Design  Support  and  Problem  Resolution  for  ASRM 

and  RSRM 

Richard  A.  Dill,  ERC  Incorporated 
R.  Harold  Whitesides,  ERC  Incorporated 

Abstract 

The  use  of  Navier-Stokes  CFD  codes  to  predict  the  internal  flow  field  environment  in  a solid 
rocket  motor  is  a very  important  analysis  element  during  the  design  phase  of  a motor  development 
program.  These  computational  flow  field  solutions  uncover  a variety  of  potential  problems  associated 
with  motor  performance  as  well  as  suggesting  solutions  to  these  problems  CFD  codes  have  also  proven  to 
be  of  great  benefit  in  explaining  problems  associated  with  operational  motors  such  as  in  the  case  of  the 
pressure  spike  problem  with  the  STS-54B  flight  motor.  This  paper  presents  results  from  analyses 
involving  both  motor  design  support  and  problem  resolution.  The  issues  discussed  include  the  fluid 
dynamic/mechanical  stress  coupling  at  field  joints  relative  to  significant  propellant  deformations,  the 
prediction  of  axial  and  radial  pressure  gradients  in  the  motor  associated  with  motor  performance  and 
propellant  mechanical  loading,  the  prediction  of  transition  of  the  internal  flow  in  the  motor  associated 
with  erosive  burning,  the  accumulation  of  slag  at  the  field  joints  and  in  the  submerged  nozzle  region, 
impingement  of  flow  on  the  nozzle  nose,  and  pressure  gradients  in  the  nozzle  region  of  the  motor. 

The  analyses  presented  in  this  paper  have  been  performed  using  a two-dimensional  axi- 
symmetric  model.  Flucnt/BFC,  a three  dimensional  Navier-Stokes  flow  field  code,  has  been  used  to  make 
the  numerical  calculations.  This  code  utilizes  a staggered  grid  formulation  along  with  the  SIMPLER 
numerical  pressure-velocity  coupling  algorithm.  Wall  functions  are  used  to  represent  the  character  of  the 
viscous  sub-layer  flow,  and  an  adjusted  k-e  turbulence  model  especially  configured  for  mass  injection 
internal  flows,  is  used  to  model  the  growth  of  turbulence  in  the  motor  port. 

The  topic  of  motor  problem  resolution  is  discussed  by  presenting  solutions  associated  with  the 
sixty-seven  second  burn  time  RSRM  motor.  The  full  motor  internal  flow  environment  for  RSRM  is 
discussed  and  the  axial  and  radial  pressure  gradients  are  shown.  The  flow  field  environment  and  pressure 
gradients  in  the  slots  are  also  discussed.  Particle  traces  from  the  burning  propellant  in  the  field  joints  are 
presented  which  show  the  tendency  of  the  center  and  aft  slots  to  collect  slag.  The  flow  field  environment 
in  the  submerged  nozzle  region  with  and  without  slag  in  the  submerged  nozzle  cavity  is  shown  and 
specific  flow  field  features  which  contribute  to  observed  post-flight  motor  erosion  patterns  is  discussed. 

The  design  support  analyses  on  the  ASRM  presented  are  for  the  zero  second  bum  time  geometry. 
The  full  motor  flow  field  environment  is  presented  along  with  axial  and  radial  pressure  gradients. 
Transition  of  the  velocity  profiles  in  the  motor  port  is  presented  and  the  effect  of  the  geometry  flare  in  the 
bore  at  the  aft  end  of  the  motor  is  shown.  The  aft  slot  deformation  analysis  is  also  presented.  This 
analysis  is  an  iterative  coupled  fluid  dynamic/mechanical  load  analysis  examining  how  two-dimensional 
flow  field  effects  in  the  motor  cause  deformation  of  the  propellant  grain.  The  submerged  nozzle  region 
flow  field  is  presented  and  discussed  as  it  relates  to  the  total  pressure  gradient  observed  in  the  aft  end  of 
the  motor.  The  radial  total  pressure  gradient  is  shown  to  be  too  great  to  allow  the  boot  cavity  motor 
pressure  measurements  to  be  compared  with  the  nozzle  end  total  pressure  computed  in  ballistic  runs. 

Conclusions  discussed  in  this  paper  consider  flow  field  effects  on  the  forward,  center,  and  aft 
propellant  grains  except  for  the  head  end  star  grain  region  of  the  forward  propellant  segment.  The  field 
joints  and  the  submerged  nozzle  are  discussed  as  well.  Conclusions  relative  to  both  the  design  evaluation 
of  the  ASRM  and  the  RSRM  scenarios  explaining  the  pressure  spikes  were  based  on  the  flow  field 
solutions  presented  in  this  paper. 
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RSRM  Motor  Geometry 


(Dashed  Line  Shows  The  67  Second  Burn  Back) 
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R5RM  FORHRRO  SLOT  AT  67  5ECQN05 
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R5RM  CENTER  SLOT  RT  67  SECONDS 
Velocity  Vectors 
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R5RM  CENTER  SLOT  RT  67  SECONDS  fluent/bf 
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RSRM  Motor  Pressure  at  67  Seconds 
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Overview 
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nozzle  assembly 
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RSRM  FULL  MOTOR  PORT  AT  0 SECONDS  BURN  TIME 
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A5RM  RFT  SLOT  AT  0 SECONDS  BURN  TIME  fluent/bf 
F i n i t e -D  i f Per  ence  Grid 20  Dc 
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TIME-ACCURATE  UNSTEADY  FLOW  SIMULATIONS  ~ 7 

SUPPORTING  THE  SRM  T+68-SEC 
PRESSURE  “SPIKE”  ANOMALY  INVESTIGATION  ^ / 

(STS-54B)  7 ' y ' 

by  A ^ 

N.  S.  Dougherty,  D.  W.  Burnette,  and  J.  B.  Holt 
Rockwell  International 
Space  Systems  Division 
Huntsville,  AL  35806 

and 

Jose  Matienzo 
MSFC/ED33 

Marshall  Space  Flight  Center,  AL  35812 
ABSTRACT 


Time-accurate  unsteady  flow  simulations  are  being  performed  supporting  the  SRM  T+68- 
sec  pressure  “spike”  anomaly  investigation.  The  anomaly  occurred  in  the  RH  SRM  during 
the  STS-54  flight  (STS-54B)  but  not  in  the  LH  SRM  (STS-54A)  causing  a momentary 
thrust  mismatch  approaching  the  allowable  limit  at  that  time  into  the  flight.  Full-motor 
internal  flow  simulations  using  the  USA-2D  axisymmetric  code  are  in  progress  for  the 

nominal  propellant  burn-back  geometry  and  flow  conditions  at  T+68-sec--Pc  = 630  psi,  y = 
1.1381,  Tc  = 6200  R,  perfect  gas  without  aluminum  particulate.  In  a cooperative  effort  with 
other  investigation  team  members,  CFD-derived  pressure  loading  on  the  NBR  and 
castable  inhibitors  was  used  iteratively  to  obtain  nominal  deformed  geometry  of  each 
inhibitor,  and  the  deformed  (bent  back)  inhibitor  geometry  was  entered  into  this  model. 
Deformed  geometry  was  computed  using  structural  finite-element  models.  A solution  for 
the  unsteady  flow  has  been  obtained  for  the  nominal  flow  conditions  (existing  prior  to  the 
occurrence  of  the  anomaly)  showing  sustained  standing  pressure  oscillations  at 
nominally  14.5  Hz  in  the  motor  IL  acoustic  mode  that  flight  and  static  test  data  confirm  to 
be  normally  present  at  this  time.  Average  mass  flow  discharged  from  the  nozzle  was 
confirmed  to  be  the  nominal  expected  (9550  Ibm/sec).  The  local  inlet  boundary  condition 
is  being  perturbed  at  the  location  of  the  presumed  reconstructed  anomaly  as  identified  by 
interior  ballistics  performance  specialist  team  members.  A time  variation  in  local  mass 
flow  is  used  to  simulate  sudden  increase  in  burning  area  due  to  localized  propellant  grain 
cracks.  The  solution  will  proceed  to  develop  a pressure  rise  (proportional  to  total  mass 
flow  rate  change  squared).  The  volume-filling  time  constant  (equivalent  to  0.5  Hz)  comes 
into  play  in  shaping  the  rise  rate  of  the  developing  pressure  “spike”  as  it  propagates  at  the 
speed  of  sound  in  both  directions  to  the  motor  head  end  and  nozzle.  The  objectives  of  the 
present  analysis  are  to:  (1)  capture  the  dynamic  responses  of  the  motor  combustion  gas 
flow  to  correlate  with  available  low-frequency  (<  12.5  sample/sec)  data  and  (2)  observe 
the  high-frequency  (up  to  50  Hz)  characteristics  of  the  response  to  determine  any 
potentials  for  dynamic  coupling. 


,i'  '--JW-’i 

paucJ837 


837 


livitwiiuNALLY  BUMK 


Rockwell 

International 

Space  Systems  Division  Huntsville  Operations 


T3 

c — 
<0  co 
~ c 
0)  o 

4"*  ■ MB 


CO  < 

c „ 


| i 

1 X 


o 

o 

X 


838 


RENG002556.03 


RockweM  STS-54B  SRM  PRESSURE 

International  "SPIKE"  ANOMALY  INVESTIGATION 

Space  Systems  Division  Huntsville  Operations 


Q 

111 

3 < o 

[[  O h 

I UJ  UJ 
03 


< 

I 

LU 


N 

N 

o 


z O 
QC  P 
LU  Z 
H Q ^ 
Z LU  UJ 
"SI 
lldh 
O 0)  t 

O)  < o 

0^0 

gs| 

3 Z 2 

2 < O 

03  : I- 

LU  O 

O * LU 
u.  ==  “3 
o ft  UJ 


UJ  = 

< OC 

£ 3 
3 CO 

O 03 

O LU 

X 

Q. 


LU 

O 

X 

o 


o 

_ H 

_ UJ 
2 03 
X z 
D O 
1.  Q. 
X 0) 
U LU 
X X 


>- 

CQ 


X 

£ 


H 

< 

—I 

LU 

QC 

QC 

O 

O 

QC 

O 


< 

H 

< 

D 

H 

X 


O 

LU 

03 


LL  <jj 

0) 

UJ 
03 
Z 

o 

CL 

a a 

o < 

2 w 

< m 

Z cm 

>-  T- 

Q uj 
uj 

QC  5 

H d 

CL  < 

< > 

O < 


03 

O 


o 0) 

< -I 


QC  < 
X Z 

O S 

^ I— 
N O 
X o. 


o >- 
m z 

O < 
I-  LU 

%i 

— QC 

o £; 

2° 

2 ° 
a h 

UJ 

£8 

x O 

O Q. 


o 


2 0. 


03  § 

in  w 


UJ 

X 


o 


LU 


LU 

> 

X 

uj  t < 

03  ^ z 

m a.  >- 

o o a 


CM 


839 


RENG002562.02 


Rockwell  STS-54B  SRM  PRESSURE 

International  "SPIKE"  ANOMALY  INVESTIGATION 

Space  Systems  Division  Huntsville  Operations 


CM 

o 

d 

to 

ir> 

C\J 

o 

o 


o 

2 

LU 

cc 


2 

§ 

1 


LU 

-I 

N 

N 

o 

z 

x 

o 

H 

o 


LU 

X 

h- 

X 

o 

LL 


o 

H 


o 

(/) 

>- 

Q 

< 

LU 

H o 

<0  LU 

1 40 
3 rs. 

lil  tf> 

X + 


CL 

o 

-I 

UJ 


(/> 

X 

o 

H 

CO 


X 

LU 

> 

o 


H H 


K- 

< 

$ 
^ o 
w -J 

Q LI- 


LLI 

CO 

I 


LU 

X 

(/> 

X 

o 

LL 

H 

< 

O 

X 


I- 

< 

CM 


o> 

CO 

X 

o 

UJ 

o 

< 

* 

o 

O 

_i 

CQ 

UJ 


O 

H 

< 

H 


O 

X 

< 


N £ 

o O 

Z «/> 

UJ  >■ 
H -J 
< < 


</)  < 


< 

UJ 

X 

< 

o 


(/> 

c/> 

Ul 


o 

H 


0) 
* 
o 
o 

X 
LU  C/> 


o 

o 


N 

N 

o 

0 z 

u X 
-J  H 

s-  ° 

J:  o 

=!  Q 

z ui 

LU  o 
Q 3 
Q Q 
3 UJ 
CO  X 

1 I 


* 

< 

UJ 

X 


< 

x 

i- 


£ 

o 


X 

UJ 

H 


X 

o 


UJ 

I- 

3 

X 

2 

O 

O 


o 

UJ 

X 

o 

< 

LU 

X 

LU 

X 

(/) 

</) 

UJ 

X 

X 


c r > 

o 

H 
C /> 

X 

UJ 

I— 

o 

< 

X 

< 

X 

o 

UJ 

> 

X 

UJ 

</> 

CQ 

0 

1 


o 

I- 

< 

X 

3 

g 

X 

z 

o 

o 

Q 

IU 

* 

o 

o 

CQ  LU 

z > 

=>  o 


UJ 

H X 

15 

UJ 
2 
o 

UJ 
O 

H UJ 

n q 


UJ 

Q 

o 

3 

0) 


N 

o 


x 

3 

h 

LU 

X 


< 

* 

o 

0 

_l 

m 

1 


>- 

x 

UJ 

> 

O 

o 

LU 

X 

o 

I- 


LU 

0) 

z 

< 

X 

I- 

£ 

o 


X 

LU 

H 


X 

o 


UJ 


O 

o 


840 


_ ,,  „ TIME-ACCURATE  CFD  ANALYSIS  OF 

international  RSRM  Pc  ANOMALY  T + 68  SEC 

Space  Systems  Division  Huntsville  Operations 


111 

H 

< 


O 

H 

X 

< 

Cl 


X 

LL 

£ 

• • 
O 

0 

o 

P 

Z 

—1 

K 

• 1 
1 

CO 

< 

_i 

_j 

< 

£ 

E 

2 

< 

a 

u. 

O 

X 

CO 

2 

K 

o 

H 

a. 

o 

co 

UJ 

§ 

§ 

co 

co 

h“ 

O 

UJ 

Li- 

CC 

i- 

< 

CD 

< 

H 

O 

UJ 

-1 

o 

2 

2 

>• 

co 

1 

UJ 

Q 

UJ 

X 

o 

a. 

< 

z 

< 

3 

co 

S 

• 

• 

• 

• 

i 

OQ 

co 

UJ 

£ 

* 

O 

o 

QC 

>- 

CQ 

Q 

UJ 

CL 

O 

-I 

UJ 

> 

UJ 

Q 

UJ 

Q 

O 

O 

< 

< O 


UJ 

h- 

< 

QC 

O 

O 

< 

■ 

UJ 


X 

UJ 

h 

QC 

O 

> 


QC 

CO 

< 


CO 

o w 

Z Q 
< D 


H 

CO 


2 

QC 
CO  z 

c o 

£ P 

o o 

u-  < 

co  DC 
co  UJ 
O H 
UJ  Z 

> 

CD  O 

UJ  H 
CO  CO 

3 3 

Z ° 

S O 
>< 

Hi 

Q 
O 
UJ 
X 
O CO 


UJ 

QC 

QC 


CO 

I- 

< 

2 

QC 

o 

UL 

Q 

UJ 

N 

a 

x 

< 

Q 

Z 

< 

h 

CO 

o 

u. 

CO 


CO 

UJ 


Q 

co 


O 

-l 

Q. 

CO 


0. 

h 


841 


RENG002355.03 


Rockwell  STS-54B  SRM  PRESSURE 

International  "SPIKE"  ANOMALY  INVESTIGATION 

Spac  Systems  Division  Huntsville  Operations 

MODELING  APPROACH 
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TEMP  3444.5  K (6200  R)  DENSITY  (GAS)  1.575E-04  LB/IN3 
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Status  of  Axisymmetic  CFD  of  an  Eleven  Inch 
Diameter  Hybrid  Rocket  Motor 

Joseph  Ruf,  ED32 
Matthew  R.  Sullivan,  EP54 
Ten  See  Wang,  ED32 
NASA/Marshall  Space  Flight  Center 
Huntsville,  AL. 


Current  status  of  a steady  state,  axisymmetric  analysis  of 
an  experimental  11"  diameter  hybrid  rocket  motor  internal 
flow  field  is  given.  The  objective  of  this  effort  is  to  develop  a 
steady  state  axisymmetric  model  of  the  11"  hybrid  rocket 
motor  which  can  be  used  as  a design  and/or  analytical  tool.  A 
test  hardware  description,  modeling  approach,  and  future  plans 
are  given.  The  analysis  was  performed  with  FDNS 
implementing  several  finite  rate  chemistry  sets.  A converged 
solution  for  a two  equation  and  five  species  set  on  a 'fine'  grid 
is  shown. 


865 


i FILMCJ 


INTENTIONALLY  BUIVK 


Status  of  Axisymmetric  CFD  Analysis  of  an  Eleven 


O 


866 


Matthew  R.  Sullivan 
Ten  See  Wang 
Marshall  Space  Flight  Center 


George  C.  Marshall  Space  Flight  Center 
Structures  and  Dynamics  Laboratory 
Computational  Fluid  Dynamics  Branch 
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George  C.  Marshall  Space  Flight  Center 
Structures  and  Dynamics  Laboratory 
Computational  Fluid  Dynamics  Branch 
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GOX  injection  pressure  = 300  psi,  flow  rate  = 6.8  Ibm/s 

GOX  injected  through  1 2 radial  ports 

overall  o/f  = 3.04 

run  duration  9.5  seconds 

nozzle  area  ratio  of  1.56 


Georgs  C.  Marshall  Space  Flight  Center 
Structures  and  Dynamics  Laboratory 
Computational  Fluid  Dynamics  Branch 
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Six  chemistry  models  tried 

- 5 species,  2 equations  up  to  1 1 species,  1 7 equations 
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VELOCITY  COLORED  BY  VELOCITY  MAGNITUDE 


ORIGINAL  PAGE  IS 

OF  POOR  QUALITY 
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UELOCITY  COLORED  BY  UELOCITY  MAGNITUDE 
Oft  Mixing  Chariber 
(ft/s) 
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N95-  23652 


Validation  of  a Computational  Fluid  Dynamics  (CFD)  Code 
for  Supersonic  Axisymmetric  Base  Flow 

P.  K.  Tucker 

NASA/Marshall  Space  Flight  Center 
Marshall  Space  Flight  Center,  AL 

Abstract 

The  ability  to  accurately  and  efficiently  calculate  the  flow  structure  in  the  base 
region  of  bodies  of  revolution  in  supersonic  flight  is  a significant  step  in  CFD  code 
validation  for  applications  ranging  from  base  heating  for  rockets  to  drag  for 
protect  ives. 

The  FDNS  code  is  used  to  compute  such  a flow  and  the  results  are  compared  to 
benchmark  quality  experimental  data.  Flowfield  calculations  are  presented  for  a 

cylindrical  afterbody  at  M = 2.46  and  angle  of  attack  a = 0.  Grid  independent 
solutions  are  compared  to  mean  velocity  profiles  in  the  separated  wake  area  and 
downstream  of  the  reattachment  point.  Additionally,  quantities  such  as  turbulent 
kinetic  energy  and  shear  layer  growth  rates  are  compared  to  the  data.  Finally,  the 
computed  base  pressures  are  compared  to  the  measured  values.  An  effort  is 
made  to  elucidate  the  role  of  turbulence  models  in  the  flowfield  predictions.  The 
level  of  turbulent  eddy  viscosity,  and  its  origin,  are  used  to  contrast  the  various 
turbulence  models  and  compare  the  results  to  the  experimental  data. 
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National  Aeronautics  and  Space  Administration 

George  C.  Marshall  Space  Flight  Center  K.  Tucker,  ED32 

Structures  and  Dynamics  Laboratory  4/21/93 

Fluid  Dynamics  Division,  CFD  Branch 


G 

o 

CD 

u 

<D 

o. 

3 

CD 

O 

MH 

0) 

O 

U 


M-l 

o 


M 

<u 


u 


3 

H 

c 

•H 

> 

a» 


>4 


880 


1993  CFD  Conference 
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Pressure-assumed  constant  static  pressure  thru  boundary  layer 

Density-calculated  via  equation  of  state 

Turbulent  kinetic  energy-interpolated  from  data  onto  grid 


INLET  VELOCITY  PROFILE 
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Third  order  upwind  Standard  k-e,  k-corr  0.034  -0.172 

Third  order  upwind  Standard  k-e,  e-corr  0.052  -0.147 

Third  order  upwind  Extended  k-e  -0.008  -0.186 

Third  order  upwind  Extended  k-e,  k-corr  0.150  -0.087 

Third  order  upwind  Extended  k-e,  e-corr  -0.179  -0.065 
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(i  + M,2)  replaces  s in  the  k-eqn  source  term 
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- Dilieate  interaction  between  compressibility  and  turbulence  generation/ transport 

- Address  vorticity  generation  issue  in  highly  compressible  flowfield 

- Complete  coarse  erid  cases  for  guidance  on  3-D  problems 
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CODE  VALIDATION  STUDY  FOR  BASE  FLOWS 

Edward  P.  Ascoli,  Adel  H.  Heiba,  Ronald  R.  Lagnado, 

Ronald  J.  Ungewitter,  and  Morgan  Williams 
Rocketdyne  Div.  /Rockwell  International 
Canoga  Park,  Ca  9 1 304 

ABSTRACT 

New  and  old  rocket  launch  concepts  recommend  the  clustering  of  motors  for 
improved  lift  capability.  The  flowfield  of  the  base  region  of  the  rocket  is  very  complex  and 
can  contain  high  temperature  plume  gases.  These  hot  gases  can  cause  catastrophic 
problems  if  not  adequately  designed  for.  To  assess  the  base  region  characteristics 
advanced  computational  fluid  dynamics  (CFD)  is  being  used.  As  a precursor  to  these 
calculations  the  CFD  code  requires  validation  on  base  flows.  The  primary  objective  of  this 
code  validation  study  was  to  establish  a high  level  of  confidence  in  predicting  base  flows 
with  the  USA  CFD  code.  USA  has  been  extensively  validated  for  fundamental  flows  and 
other  applications.  However,  base  heating  flows  have  a number  of  unique  characteristics 
so  it  was  necessary  to  extend  the  existing  validation  for  this  class  of  problems. 

In  preparation  for  the  planned  NLS  1 .5  Stage  base  heating  analysis,  six  case  sets 
were  studied  to  extend  the  USA  code  validation  data  base.  This  presentation  gives  a 
cursive  review  of  three  of  these  cases.  The  cases  presented  include  a 2D  axi-symmetric 
study,  a 3D  real  nozzle  study,  and  a 3D  multi-species  study.  The  results  of  all  the  studies 
show  good  general  agreement  with  data  with  no  adjustments  to  the  base  numerical 
algorithms  or  physical  models  in  the  code.  The  study  proved  the  capability  of  the  USA 
code  for  modeling  base  flows  within  the  accuracy  of  available  data. 
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MORGAN  WILLIAMS 

ROCKETDYNE  DIV/  ROCKWELL  INTERNATIONAL 
MFSC  CFD  WORKSHOP  APRIL  20-22  1993 


BASE  HEATING  STUDY  VALIDATION 
OBJECTIVE  AND  APPROACH 
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Rockwell  International 
Rocketdyne  Division 


SIX  CASE  SETS  COMPLETED  TO  EXTEND  USA  CODE 
VALIDATION  FOR  BASE  HEATING  ANALYSIS 
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Rockwell  International 

Rocketdyne  Division 


AGARD  NOZZLE  AFTERBODY  CFD  VALIDATION  STUDY 
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AGARD  CFD  MODELING  APPROACH 
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PERFECT  GAS 

BALDWIN-LOMAX  TURBULENCE  MODEL 


AGARD  CFD  VALIDATION  STUDY  CONCLUSIONS 
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OVERALL  AGREEMENT  APPEARS  TO  BE 
GOOD  AND  THE  KEY  FLOW  FEATURES 
AT  THE  BASE  ARE  CAPTURED 


AEROSPIKE  CFD  VALIDATION  STUDY 
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Kingsland,  R.B.,  Petrilla,  S.,  And  Baker,  W.M.,  "SSTO  Aerospike 
Nozzle  Demonstrator  Test  in  the  Fluidyne  Channel  No.  9 
Test  Facility,"  SSTO  Pretest  Conference,  March  1991. 


AEROSPIKE  CFD  MODELING  APPROACH 
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BALDWIN-LOMAX  TURBULENCE  MODEL 


AEROSPIKE  CFD  ANALYSIS  RESULTS 
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AEROSPIKE  CFD  ANALYSIS  RESULTS  (CONT’D) 

LOW  PRESSURE  RATIO  CASE 
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UVA  CFD  VALIDATION  STUDY 
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J.  McDaniel,  D.  Fletcher,  R.  Hartfield  and  S.  Hallo,  "Staged  Transverse  Injection 
Into  Mach  2 Flow  Behind  a Rearward-Facing  Step:  A 3-D  Compressible  Test 
Case  for  Hypersonic  Combustor  Code  Validation,"  AIAA-91-5071 . 
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• O-EQUATION  GOLDBERG  TURBULENCE  MODEL 


UVA  CFD  ANALYSIS  RESULTS 
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CFD  CODE  VALIDATION  ASSESSMENT 
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FOUR-NOZZLE  BENCHMARK  WIND  TUNNEL 
MODEL  USA  CODE  SOLUTIONS  FOR  SIMULATION  OF  MULTIPLE 
ROCKET  BASE  FLOW  RECIRCULATION  AT  145,000  FT  ALTITUDE 

by 

N.  S.  Dougherty  and  S.  L.  Johnson 

Rockwell  International 
Space  Systems  Division 
Huntsville,  AL  35806 


ABSTRACT 


Multiple  rocket  exhaust  plume  interactions  at  high  altitudes  can  produce  base  flow 
recirculation  with  attendant  alteration  of  the  base  pressure  coefficient  and  increased  base 
heating.  A search  for  a good  wind  tunnel  benchmark  problem  to  check  grid  clustering 
technique  and  turbulence  modeling  turned  up  the  experiment  done  at  AEDC  in  1961  by 
Goethert  and  Matz  on  a 4.25-in.  diameter  domed  missile  base  model  with  four  rocket 
nozzles.  This  wind  tunnel  model  with  varied  external  bleed  air  flow  for  the  base  flow  wake 
produced  measured  p/pref  at  the  center  of  the  base  as  high  as  3.3  due  to  plume  flow 
recirculation  back  onto  the  base.  At  that  time  in  1961,  relatively  inexpensive 
experimentation  with  air  at  y = 1.4  and  nozzle  Aq/A*  of  10.6  and  0n  = 7.55  deg  with  Pc  = 
155  psia  simulated  a LO2/LH2  rocket  exhaust  plume  with  y=  1.20,  Ae/A*  of  78  and  Pc 
about  1,000  psia.  An  array  of  base  pressure  taps  on  the  aft  dome  gave  a clear 
measurement  of  the  plume  recirculation  effects  at  poo  = 4.76  psfa  corresponding  to 

145,000  ft  altitude.  Our  CFD  computations  of  the  flow  field  with  direct  comparison  of 
computed-versus-measured  base  pressure  distribution  (across  the  dome)  provide 
detailed  information  on  velocities  and  particle  traces  as  well  eddy  viscosity  in  the  base 
and  nozzle  region.  The  solution  was  obtained  using  a six-zone  mesh  with  284,000  grid 
points  for  one  quadrant  taking  advantage  of  symmetry.  Results  are  compared  using  a 
zero-equation  algebraic  and  a one-equation  pointwise  Rt  turbulence  model  (work  in 
progress).  Good  agreement  with  the  experimental  pressure  data  was  obtained  with  both; 
and  this  benchmark  showed  the  importance  of:  (1)  proper  grid  clustering  and  (2)  proper 
choice  of  turbulence  modeling  for  rocket  plume  problems/recirculation  at  high  altitude. 
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RnfWli  FOUR-NOZZLE  CLUSTER  HIGH  ALTITUDE 

International  BASE  FLOW  BENCHMARK 

Space  Systems  Division  Huntsville  Operations 
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NO-GRADIENT  OUTLET 

BOUNDARY  OUTER  RADIUS  IS  A NO-SLIP  WALL 

BOUNDARY  (OMITTED  FROM  FIGURE  FOR  CLARITY) 


MRCH  NUMBER 
Four  Nozzle  Base  Flow 
Altitude  =145  kft,  Pref  =4.76  psfa 
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932  ORIGINAL  PAGE  IS 

OF  POOR  QUALITY 


FOUR  NOZZLE  BASE  PRESSURE 
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Base  Radius  (in) 
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PARTICLE  TRACES 
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Ui'ilGinAL  ^AGE  IS 
OF  POOR  QUALITY 


UELOC I TV  COLORED  BV  flflCH  NUMBER 
Four  Nozzle  Base  Flow 
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NUMERICAL  STUDY  OF  BASE  PRESSURE  CHARACTERISTIC  CURVE 
FOR  A FOUR-ENGINE  CLUSTERED  NOZZLE  CONFIGURATION 

Ten-See  Wang 

Computational  Fluid  Dynamics  Branch 
NASA  - Marshall  Space  Flight  Center 
Marshall  Space  Flight  Center,  AL  35812 

Abstract 


Excessive  base  heating  has  been  a problem  for  many  launch 
vehicles.  For  certain  design  such  as  the  direct  dump  of  turbine 
exhaust  in  the  nozzle  section  and  at  the  nozzle  lip  of  the  Space 
Transportation  Systems  Engine  ( STME ) , the  potential  burning  of  the 
turbine  exhaust  in  the  base  region  have  caused  a tremendous 
concern.  Two  conventional  approaches  have  been  considered  for 
predicting  the  base  environment:  (1)  empirical  approach,  and  (2) 
experimental  approach.  The  empirical  approach  uses  a combination 
of  data  correlations  and  semi-theoretical  calculations.  It  works 
best  for  linear  problems,  simple  physics  and  geometry.  However, 
it  is  highly  suspecious  when  complex  geometry  and  flow  physics  are 
involved,  especially  when  the  subject  is  out  of  historical 
database.  The  experimental  approach  is  often  used  to  establish 
database  for  engineering  analysis.  However,  it  is  qualitative  at 
best  for  base  flow  problems.  Other  criticisms  include  the 
inability  to  simulate  forebody  boundary  layer  correctly,  the 
interference  effect  from  tunnel  wall,  and  the  inability  to  scale 
all  pertinent  parameters.  Furthermore,  there  is  a contention  that 
the  information  extrapolated  from  subscale  tests  with  combustion 
is  unconservative. 

One  potential  alternative  to  the  conventional  methods  is  the 
computational  fluid  dynamics  (CFD) , which  has  none  of  the  above 
restrictions  and  is  becoming  more  feasible  due  to  maturing 
algorithms  and  advancing  computer  technology.  It  provides  more 
details  of  the  flowfield  and  is  only  limited  by  the  computer 
resources.  However,  it  has  its  share  of  criticism  as  a predictive 
tool  for  base  environment.  One  major  concern  is  that  CFD  has  not 
been  extensively  tested  for  base  flow  problems.  It  is  therefore 
imperative  that  CFD  be  assessed  and  benchmarked  satisfactorily  for 
base  flows. 

In  this  study,  the  turbulent  base  flowfield  of  a experimental 
investigation  for  a four-engine  clustered  nozzle  is  numerically 
benchmarked  using  a pressure  based  CFD  method.  Since  the  cold  air 
was  the  medium,  accurate  prediction  of  the  base  pressure 
distributions  at  high  altitudes  is  the  primary  goal.  Other  factors 
which  may  influence  the  numerical  results  such  as  the  effects  of 
grid  density,  turbulence  model,  differencing  scheme,  and  boundary 
conditions  are  also  being  addressed.  Preliminary  result  of  the 
computed  base  pressure  agreed  reasonably  well  with  that  of  the 
measurement.  Basic  base  flow  features  such  as  the  reverse  jet, 
wall  jet,  recompression  shock,  and  static  pressure  field  in  plane 
of  impingement  have  been  captured. 


Numerical  Study  of  Base  Pressure  Characteristic  Curve 
for  a Four-Engine  Clustered  Nozzle  Configuration 
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Base  environment  predictive  methods 


© 

E © 
o o 

© © 
O)  > 

4.  * 

e 3 


(0 

TJ 

C 

© 

ws 

(0 

o 


(0 

W 

E 

0) 

.O 

O 


CO 

0) 

c 


o 

CO 

o 

u 

Q.  — 

8-  o 

8 8 

V-  _Q 

'o.  0) 
E ■* 

© O 

£ 5 


CL 

E 

o 

o 

■o 

c 


© « 
> £* 
Q.  © 

© E 

Q.  ° 

E 


© 

O) 


© 

Q. 

E 

o 

o 

c 

© 


© 

3 

o 

■ mam 

o 

a mrnm 

Q. 

© 

3 

© 

O) 


TJ 

© 

> 

O 

> 

c 

a^M 

© 

i_ 

© 

© 

O 

© 

© 

© 

.Q 

© 

© 

O 

3 

© 


© 

© 

© 

A 

© 

© 

■o 

© 

o 


© 


3 

O 

© 
a as 

O 

© 

A 

3 

© 

© 


C 

© 


© 

o 

© 

CL 

© 

© 


© 

c 

© 

cn 

c 

© 

© 

c 

'5> 

c 

© 


o 

© 

O 

Q. 

a 

© 

© 

c 

© 

E 

© 

a 

x 

© 

© 


© 

© 

© 

A 

© 

© 

■o 

© 


© 

BMB 

A 

© 

© 

© 


TJ 

© 

© 

3 

C 

© 

o 


© 

© 

© 

c 

> 

To 

o 

© 

£ 

a^M 

*-> 

© 

© 

O 

a mb 

© 

c 

Q. 

Q. 

TJ 

C 

c 

3 

4-* 

© 

O 

E 

£ 

.Q 

o 

O 

> 

H— 

t- 

■o 

© 

o 

o 

© 

© 

© 

© 

9= 

A 

o 

© 

u. 

© 

O 

© 

o 

c 

© 

© 

© 

© 

3 

© 

.Q 

E 

t: 

4-* 

© 

w 

© 

© 

o 

c 

8MB 

> 

BUB 

4-* 

> 

© 

© 

♦4 

13 

MB 

aBBB 

a^M 

© 

a mb 

A 

gmt 

© 

© 

3 

© 

o 

O' 

L. 

■ M 

a 

i 

1 

■ 

* 


944 


- inability  to  scale  all  pertinent  parameters 

- information  extrapolated  from  subscale  test  with  combustion  is 
unconservative 
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0 Two-Equation  Turbulence  Model 


Parametric  Study 
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distance  from  center  of  heat  shield,  r.  Inches 
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★ The  potential  of  using  CFD  as  a predictive  tool  for 
base  environment  prediction  is  demonstrated 
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